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Abstract 

Using extended Khuri-Treiman equations, we evaluate the final state in- 
teractions due to two-pion rescatterings to the decays r\ — ► 7r°7r + 7r _ and 
r\ —* 7r 7r 7r°. As subtraction to the dispersion relation we take the one- loop 
chiral perturbation theory result of Gasser and Leutwyler. The calculated 
corrections are moderate and amount to about 14% in the amplitude at the 
center of the decay region. A careful analysis of the errors inherent to our 
approach is given. As a consequence, the experimental rate of the decay can 
only be reproduced if the double quark mass ratio Q~ 2 = " ld ~ m :" • is 
increased from the usual value of 1/(24. 1) 2 to 1/(22.4 ± 0.9) 2 . We have also 
calculated the ratio of the rates of the two decays and various Dalitz Plot 
parameters. In particular, the linear slope a in the charged decay is different 
from the one- loop value and agrees better with experiment. 
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1 Introduction 



Chiral perturbation theory (ChPT) [[T], [2| offers a consistent description of low 
energy QCD, in particular of the strong and weak mesonic interactions. Whereas 
most phenomena are quite well accounted for, the predictions for the decays 77 — > 3n 
remain much below the experimental results. The lowest order amplitude is fully 
determined by chiral symmetry || to be 

= -M^Zpl. (1.1) 

/< 2 >( s „) = rw = i + %4 

m~ — 7774 

where m u and are the up and down quark masses, F and B QCD parameters 
and s a , Sb, s c , s kinematical variables to be defined later. While the form of eq. ( |1.1| ) 
is valid for the decay into 7r°7r + 7r~, the expression for three neutral pions is obtained 
by setting f^ 2 '(s a ) = 1. Furthermore, in lowest order in the chiral expansion one has 

B (m d -m u ) = (m 2 K0 - m 2 K+ ) - (mlo - ml + ) (1.2) 
= ml. 

These equations yield a width of 66 eV for the decay 77 — > 7r°7r + 7r~ which is far below 
the experimental value of 281 ± 28 eV [|J. It is an unusual situation, since most 
predictions of current algebra fall within 20 to 30% of the experimental results. 

However, also the next-to-leading calculation of Gasser and Leutwyler || failed to 
reproduce the experimental value, despite a dramatic improvement. These authors 
obtain a width of 167 ± 50eV for the 77 — > 7r°7r + 7r~ decay (We have adopted the 
'new' value of f w |J which increases the width by about 7 eV). Thus, one may con- 
clude that higher order corrections are very large, possibly preventing a satisfactory 
representation of this decay by ChPT. 

Since the decay rate is proportional to (m u — m^) 2 it is particularly sensitive to the 
value of the quark masses. The quark mass ratios are known with some precision [^] 
from a variety of low energy investigations. Since they are fundamental parameters 
of the basic theory, they should be determined as accurately as possible. Of special 
interest is the up quark mass because m u = is a very appealing solution to the 



strong CP problem [g] . From eq. ( |1 . 1| ) we see that decreasing m u indeed increases 
the rate as required by experiment. On the other hand, corrections to the second 
factor will also have this effect. Thus, only a careful investigation of higher order 
effects will enable one to draw a conclusion on the quark masses. 

Let us formulate this more precisely. The amplitude for the decay 77 — > 7r°7r + 7r~ 
can be written as 

Bo(md — m u ) 

F(s a , s b , s c ) = 3 ^ ?2 f{s a , s b , s c , A, m u , ..) + e.m. (1.3) 
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where the quark masses are the renormalization group invariant masses, A a QCD 
scale, s a , Sb, s c the Mandelstam variables for the decay. Sutherland's theorem [|] 
implies that the electromagnetic contribution is of order p 2 , where p 2 stands for any 
invariant product of momenta. Moreover, the leading (in momentum) electromag- 
netic contribution of order is further suppressed by a factor if the amplitude 

is assumed to be linear in s a [ |TC| , |TTf . Recently, these terms were reanalyzed in an 
effective Lagrangian framework []I2] and the expected smallness of the electromag- 
netic corrections (Sutherland's theorem) was confirmed. We will therefore neglect 
them throughout. 

Returning then to the QCD contribution in eq. ( |1.3|) , the second factor, the 
function /, is expanded in chiral perturbation theory in powers of momentum and 

mass 

/ = / (2) + / {4) + (1-4) 

Here, has already been defined in eq. ( |1 . 1| ) and etc. are of higher orders. 

was given explicitely by Gasser and Leutwyler || and consists of a variety of 
loops and counterterms. The corresponding coupling constants also occur in other 
observables which are also calculated using the chiral Lagrangian, for instance the 
decay constants or the masses of the mesons. To order p 4 , it is therefore possible 
to express some of these corrections through physical quantities. In particular, the 
quark mass contribution to the kaon mass difference can be written as |[T3| 
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AK QCD = (m 2 K0 - m 2 K+ ) QCD = Q- 2 ^( m 2 K - ml) (1.5) 

and can be obtained from the measured masses once the electromagnetic mass dif- 
ferences are subtracted. Here, 

2 m d -m u m d + m u 

Q = ; — - (1-6) 

m s — m m s + m 



with m = \{m<i + m u ). 

Then, one can rewrite the amplitude as 



/•K.-s/,.*.) -Q- 2 ^(m 2 K -ml)^-^ (1.7) 
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■f (2) (s a )(l + 5 DEC ). 

(5 dec contains the remaining part of the 0(p A ) terms as well as corrections from 
higher orders in the chiral expansion. Terms of order (m^ — m u ) 2 are tiny and are 
therefore neglected throughout. As anticipated, the rate is indeed proportional to 
the quark mass difference squared. Q 



In passing, we note that in the generalized framework of chiral perturbation theory 1 14 



are more complicated and a careful reanalysis of all relations entering eq. (1.7) is required. 



matters 
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The decay 77 — ^ 37T fixes the quantity Q 2 rather than the ratio This is 



in agreement with the observation by Kaplan and Manohar Jl5| that the chiral 



Lagrangian possesses a curious reparametrization invariance under which masses 
and some of the other parameters are changed. It implies that matrix elements 
derived from Green functions which do not involve scalar or pseudoscalar densities 
are invariant |1(J with respect to that transformation. And indeed, the combination 
Q 2 in eq. ( |1.6| ) is invariant, while for instance ^ is not. Unfortunately, there is 
no direct experimental access to this ratio, the relevant Green functions must be 
obtained from semiphenomenological analyses [0]. 

There exists an independent determination of Q 2 by means of the following equa- 
tion 

ml = AKqcd + (AK - A7r) em . (1.8) 

Provided that the electromagnetic contribution is known, eqs. ( |1.5| ) and ( |1.8j ) can 
be used to check the consistency of the different calculations or to predict the rate 
for the decay r\ — > 3n. 

The calculation of the electromagnetic mass differences is notoriously difficult. 
Dashen's theorem |T7| implies that the electromagnetic mass difference of the kaons 
is equal to the mass difference of the pions, [] up to terms of order e 2 p 2 . If these are 
neglected, one has AKqcd = ml and one can write 

F(s a , s b , s c ) = ~^r 2 ■ f (2) (s a )(l + 5 D ec) (1.9) 

and 5 dec determines the rate. In their calculation of the decay rate for 77 — ► 3n 
Gasser and Leutwyler || obtained a correction 5dec of about 50%. Together with 
eq. ( |1.9|) , the value of 167 eV for the width resulted, still far below the experimental 
number. 

However, it was found recently [18|, [19[] that the corrections to Dashen's theorem 



could be substantial. One should therefore start from eq. (|1.7|), avoiding unnecessary 
uncertainties arising from corrections to Dashen's theorem. The normalization of 
the 77 — > 37r amplitude is then determined by Q 2 which follows in principle from eq. 



( |1 . 5| ) . As pointed out in pD |, the positive corrections to Dashen's theorem found in 
JL8| increase AK such that together with the value of 5 dec obtained by Gasser and 
Leutwyler the experimental rate for 77 — > 3tt would be reproduced. 

Thus, the rate of the 77 decay can be increased either by determining further cor- 
rections in 5 dec (beyond those calculated by Gasser and Leutwyler) or by increasing 
the value of AKqcd- in view of eq. (|1.6| ) this latter term amounts to an increase in 
Q 2 and thus to a smaller up quark mass. But before any conclusion on the quark 
masses can be drawn, the decay corrections must be understood better. 

It has been argued [^TJ that 777/ mixing may increase the theoretical rate to the 



experimental value. However, it was recently shown [22] that the effects of the rj' 



are in fact included in the standard treatment and that the enhancement is due to 



Hhe QCD contribution to the pion mass difference is proportional to {rrid — m u ) 2 and negligible 
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an incomplete treatment of the resonances as an explanation for the constants of 
the chiral Lagrangian. We therefore omit this issue in the following. 

The goal of this paper is to calculate one class of corrections in 6dec, the so called 
unitary corrections which ensure that the decay amplitude satisfies unitarity. These 
corrections are the sum of certain diagrams, namely those describing the rescattering 
of the final state particles to all orders. It is believed that they dominate the 
corrections, in particular since the final state interactions of two pions in the 1 = 
s-wave channel is strong and attractive. In fact, in the one- loop calculation of Gasser 
and Leutwyler they yield the largest new contributions. They account for roughly 
three quarters of the corrections, pion loops contributing about 85% and kaons and 
etas about 15% thereof, while other terms provide the remaining quarter (of course, 
the division of the corrections depends on the renormalization scale chosen; however, 
physical arguments such as resonance saturation of the counterterms indicate that 
the scale used for the above numbers is reasonable). 

A method to evaluate the unitary corrections has been given long ago by Khuri 
and Treiman |23J]. Keeping only two particle rescatterings of pions (and neglecting 
p-waves) they derived a set of dispersion relations, the Khuri- Treiman equation. 



This paper was subsequently analyzed by several authors [ER Eq, E7L E8L 29, pQ]. In 



particular, Kacser [28| derived the correct prescription for the analytical continuation 
of the partial wave projections of the decay amplitude. Bronzan |29| and later also 
Neveu and Scherk omitting also the 1 = 2 rescatterings, solved the Khuri- 
Treiman equation approximately. Using the properties of the Omnes function D(s) 
entering their approximate solution they argued that the remaining terms are small. 

Some time ago, Roiesnel and Truong |32] reconsidered the problem, using similar, 
although not equivalent methods to those of Neveu and Scherk to calculate the 
unitary corrections. These authors found a large enhancement of the amplitude, 
enough to account for the observed rate. The numerical results of this work were 
criticized in Ref. |J on physical grounds and we confirm their findings. We believe 
therefore that their result is a overestimate of the unitary corrections. 

In this work we will use (and generalize by the inclusion of p- waves (see also 
the method of Khuri and Treiman to determine numerically the unitary corrections. 
It is organized as follows. 

In section 2 we review the general structure of the decay amplitude for 77 — > 37r and 
give the necessary definitions. In section 3 we discuss the Khuri- Treiman equations 
and the inclusion of the 1 = 1 rescatterings. In section 4 we discuss the subtraction 
procedure, and in section 5 we build up an iterative scheme for the numerical solution 
of the set of coupled equations for the projected amplitudes U, V, W. In section 6 
we discuss our numerical results and the quality of our iteration, and then in section 
7 we draw the conclusions. Some appendices contain necessary technical details. 
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2 Isospin decomposition 



In this section we discuss the general form of the amplitude for the decay rj — > 3tt. 
It serves to fix the notation, to give the isospin decomposition of the amplitude and 
to discuss its symmetry properties. 

We denote the three pions in the decay rj — > 3ir by the letters a,b,c, their cor- 
responding four-momenta by k a , kb, k c and label their isospins by a, (3, 7. The 
invariant masses of the pions are then k 2 a = k\ = k 2 c = ml disregarding all mass 
differences at this stage. The four momentum of the isospin singlet rj is denoted by 

r- 

In the standard model of the strong, electromagnetic and weak interactions, the 
decay 77 — ^ 37T proceeds either through the isospin breaking piece of the QCD hamil- 
tonian or through an operator of electromagnetic origin. For both contributions, 
the interaction Hamiltonian H(0) is AI = 1. It is then convenient to treat the 77 
as a spurious isospin triplet and to assume isospin conservation in the decay. The 
corresponding isospin state is labeled by p. 

The amplitude for a general isospin assignment of the pions is 

M = < a,k a ;P,k b ;-f,k c \H(0)\k r > 

= i(2ir) 4 5(k a + k b + k c - £v)M Q/37;p (s a , s b , s c ). (2.1) 

The invariant amplitude has the decomposition 

M a p T:P (s a , s b , s c ) = F(s a , s b , s c )5^5 ap 

+F(s b , s c , s a )8 ai 8p p 

+F(s c ,s a ,s b )5 a(3 5 7P , (2.2) 

where the isoscalar amplitude F is a function of the Mandelstam variables 

Sa = {k r ka) , Sb — [kr kb) , s c = (k r k c ) . (2-3) 

Of course only two of these variables are independent as they are related by the 
on-shell condition 

s a + s b + s c = 3ml + m^ = 3s . (2.4) 
Bose statistics of the three pion system implies the symmetry 

F(s a ,s b ,s c ) = F(s a ,s c ,s b ) (2.5) 

such that the amplitude M a p r ^ p (s a , s b , s c ) remains invariant under the exchanges of 



The amplitude for the decay into charged pions is then found to be 

< 7T+, 7T- 7r°|#(0)|77 >= t(2n) 4 5(k a + k b + k c - k r )F(s a , s 6 , s c ), (2.7) 

whereas the amplitude for the decay into neutral pions is expressed as 

< 7T°, 7T°, 7r°|iy(0)|?7 > = i(27i) 4 5(k a + k b + k c - k r ) (2.8) 

[F(s a , s b , s c ) + F(s b , s c , s ) + F(s c , s , s 6 )] . 

3 Isospin sum and scalar dispersion relation 

We begin with the dispersion representation for the decay amplitude eq. (|2 
derived by Khuri and Treiman 



7T J4 S' — S„ — IE 



i /-oo disci, Mp ja .Js a , s' b ,s c ) 
+ ~ ds' b ^ ; (3.1) 

7T J4 S fc — Sft — IE 

+ I Z" 00 disc Ic M 7»/3;p( g q ; g ^ 4) 



7T J4 S' c — S c — IE 

Eqn. (|3.1|) expresses the connection between causality and analyticity for the consid- 
ered amplitude. Unitarity gives an expression for the discontinuity of M a/3rtP in s a 
for fixed s c and analogous ones for the other two discontinuities in the form of a sum 
over intermediate states involved in all the rescattering processes. The intermediate 
state with lowest mass contributing in this sum is the two-pion state leading to a 
rescattering of two outgoing pions in the decay. As any other intermediate state 
such as a KK or a Air state will contribute to the discontinuity at a much larger 
threshold it is a reasonable approximation to drop them |24j|. Hence, we find 



2% disc | a M Q/37 . p (s a , s 6 , s c ) = J d Pddp e (3.2) 

' 8,£ 

(27r) A 5(k r -k a -p d - p e )N* S£ . Pl {s a , s' d , s'J • M a5£ . p (s a , s d , s e ) 

where p a and p e are the four-momenta of the intermediate pions d, e. In eq. (|3.2| ) 
we have introduced the invariant pion-pion scattering matrix element 

< a,k a ;(3,k b \T\E,p e ;r],pf > (3.3) 
= i(2Tc) 4 5(k a + k b -p e - p f )N a/3;£V (s a , s' d , s'J 

and the phase space measure for the intermediate pions 
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The variables s' d and s' e are defined by 

s' d = {h-Pdf (3.5) 



s'e = {h-Pt 



\2 



Note also that all momenta are taken in units of m n . In the Khuri-Treiman approx- 
imation, the discontinuity of the decay amplitude becomes thus just a product of 
the pion scattering amplitude with the decay amplitude itself. 

To perform the phase space integrals in eq. (|3.2|), one usually chooses the c. m. 
frame of the rescattered pions defined by + = 0. Fixing the #-angle in dpd to 
be the angle 9d a between and p d the discontinuity becomes 



disc | a M Q/37;p = _L/3( Sa ) (3.6) 

OZ7T 

da 



Yl J -^-N&;/*r( s a,COS0 M ) • M aBe .p(s a ,CO&0 dl ?, 



S,e 

where we have introduced the shorthand notation 



M = (3-7) 

We note that some of the kinematical variables are not well behaved in the decay 
region; we shall supplement a correct analytic interpretation for them below. 

Adopting the normalizations from Ref . |54 , we continue with the isospin decom- 
position of the invariant pion-pion scattering matrix element 

Nafc-ysiSa, Sb, s c ) = ~(T°(s a , s b , s c ) - T 2 (s a , s b , s c ))5 a/3 5 lS 

+ ^(T 2 (s a ,s b ,s c ) + T\s a ,s b , s c ))6 ay 6ps (3.8) 
+ 7^(T 2 (s a , s b , s c ) - T x (s a , s b , s c ))5 a5 5 M 

where the T 1 are the isospin amplitudes to isospin /. In our frame of reference 
they depend on two variables only, the total invariant mass of the pions and the 
intermediate scattering angle 9 b d- The corresponding isospin decomposition of the 
decay amplitude has already been given in eq. (|2.2|) . 

We now insert eq. ( |3.8| ) into eq. ( |3.6| ) and use eq. ( |2.2| ) This yields for the 
discontinuity of M in s a 



1 „, , rdtt 



disc | a M^ 7;p = 3^:^(sa) J 



da 



An 

[Sf3yS ap {T°*(s a , s' d , s' e )F(s a , s d , s e ) 
+ i(T°*( Sa , s' d , s' e ) - T 2 *(s a , s' d , s' e ))F(s e , Ba , s d ) 
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0*f„ I >\ 



T *(s a ,4,s' e ))F(s d ,s e ,s a )} 



1 



+5 7a 5 /3p {^(T 2 *(s a , s^, s'J + T u (s a , s' d , s' e ))F(s d , s e , s a ) 
(T 2 *(s a , s' d , s' e ) - T u (s a , s' d , s' e ))F(s e , s a , s d )} 



(3.9) 



+5 aP 5 lp {-{T 2 *(s a) s' d , s' e ) + T 1 *(s a , s' d , s' e ))F(s e , s a , s d ) 

+ i(T 2 *(s a , s' d , s' e ) - T 1 *^, s' d , s' e ))F{s d , s e , s a )}}. 

Next, we use this result and the two analogous ones for the discontinuities in s b 
and s c in the dispersion representation eq. (|3.1|) for the full amplitude. Comparing 
the coefficients in the different isospin channels we obtain three dispersion relations 
for the scalar functions F [ f24|j : 



F{s a , Sb, s c 



ds' 



7T J 4 S' a — >. , , — %E 327T 
-i0* I J J J 



dtt 



da 



47T 



{T *{s' a ,s' d ,s' e )F(s' a ,s d ,s e ) 
+ l(T°*(C s'd, O - T 2 *«, s' d , s' e ))F(s d , s e , s' a ) 

+l(T°*(s' a , s' d , s' e ) - T 2 *«, s' d , s' e ))F(s e , s' a , Sd)} 



ds' h 



7T J 4 S, 



Sb — is 327T 



dtt 



db 



An 



(3.10) 



{-(T 2 *(4, s' d , s' e ) + T 1 *^, s' d , s' e ))F(s e , s' b , s d ) 



1 



(T 2 *(s' b , s' d , s' e ) - T u (s' b , s' d , s' e ))F(s d , s e , s' b )} 



+ - 



ds' 



7T J 4 S'—S r — IE 327T 



d^dc 
An 



{^(T 2 *(s' d , s' e , s' c ) + T 1 ^, s' d , s' e ))F(s d , s e , s' c ) 

+ \(T 2 *(s' c , s' d , s' e ) - T 1 *( S ' C) s' d , s' e ))F(s e , s' c , s d )}. 

The invariant variables s d , s e , s' d , s' e must be expressed in terms of s a , Sb, s c . The other 
two scalar dispersion representations are obtained from eq. (|3.10|) by permutation 
of a, b, c. 

To evaluate further this dispersion relation we retain only the s- and p-wave 
contributions to the rescattering. The relevant partial wave expansions of the pion 
scattering amplitudes T 1 for fixed isospin / are given as usual by 



32-7T 00 

I^cosfl) = ^£(2Z + l)P,(cos0)if( S ) 



(3.11) 
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where s is the total invariant mass and 9 the scattering angle in the cm. system of 
the colliding pions. In the elastic region 4 < s < 16 the partial wave amplitudes are 
parametrized by real phase shifts as 



/f( S ) = e^>sin#( S ). 



The lowest contributions for 7 = 0,2 are s-waves 



T 2 (s) 

while for I = 1 it is a p-wave 

T 1 (s,cos9) 



32tt 
W) 



(3.12) 



(3.13) 



32tt 



fo(s) 



32tt 
W) 



3cos# • fl(s). 



(3.14) 



The factors -^(3 in ( |3.10|) just cancel now against the corresponding factors in the 
partial wave expansion of T 1 and we are left with two types of angle integrals. 

In the case of s-wave contributions the pion-pion rescattering is isotropic. For 
I = 0, 2 we have thus angle integrations of the type p3 



F(s a ) = 

F(Sa) = 



da 



du 



47T 
dQda 

Air 



F(s a ,s d ,s e ) 
F(s d ,s e ,s a ) 



(3.15) 



where Sd, s e are expressed through s a and cos 9^ as 



Sd 
d^da 



3s - s a - s d 
3s — s 



- + K(s a ) ■ cos9 da 



(3.16) 



dip ■ d cos 9 da . 



Here we introduce the Kacser function 2< 



K(s) 



[s - 4)(s - (m - l) 2 )(s - (m + l) 2 ) 



4s 



(3.17) 



Its cuts and limits in the complex s-plane will be discussed below. 

We may now recast the integrations above as integrations over Sd and obtain 



F(s a ) 

and in an analogous way 



2K(s a ) Js d „( Sa ) 



(Sa) 



ds d F(s a , s d , 3s - s a - s d ) 



F(s r , 



2K(s a ) J Sd _(s a ) 



Sd+(Sa) 



ds d F(s d , 3s - s a - s d , s a ) 



(3.18) 



(3.19) 
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where 

s d+ (s)^^^ + K(s), 8 ^ 8 ) = ^1-K(8). (3.20) 
For later use we note a symmetry property of the integral F(s a ) 

F ( s a) = n7 w \ / ds d F(3s - s a - s d , s a , s d ) (3.21) 

ZK {S a ) Js d -(s a ) 

which follows immediately from the Bose symmetry properties of F itself and a 
change of variable s d — > 3s — s a — s d . 

The p-wave contribution for / = 1 leads to the more complicated integral 

I(s a ,s b ,s c ) = J cos 9 bd F(s e , s a ,s d ). (3.22) 

In order to perform the angle integration explicitly and to get rid of the two angles 
we invoke a partial wave expansion for F 

oo 

F(s e ,s a ,s d ) =^(2/ + l)P,(cos^ a )^ d (sa). (3.23) 



Making use of well-known properties of the Legendre polynomials we obtain 

I(s a , s b , s c ) = cos6 ba ■ F? d (s a ). (3.24) 
The I = 1 projection of F is next expressed as an integral over s d 



1 rSd+(Sa) 

AK 2 (s a ) Js d _(s a ) 

■F(3s - s a - s d , s a , s d ) 



F i d ( s «) = JTT^ r +(Sa) ds d {2s d + s a -3s ) (3.25) 

4A (Sa) Js d _(s a ) 



and cos 6 ba is recast as 



cos9ba = mty (3 - 26) 



Plugging all together we obtain the desired expression 



I(s a , s b , s c ) = ^ S \ f + ds d (2s d + s - 3s ) (3.27) 
8K'\s a ) Js d _(s a ) 



•F(3s Q - s a - s d , s a , s d ). 
Let us note for later use a symmetry property of I(s a , s b , s c ) 



I{s a ,s b ,s c ) = , S< \ / ds d (2s d + s a - 3s ) (3.28) 

oK d (s a ) Js d _(s a ) 

■F(s d ,3s - s a - s d ,s a ) 
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following from the Bose symmetry properties of F and a change of variable Sd — > 
3so — s a — sa- Finally we decompose I(s a , s&, s c ) into 

I(s a ,s b ,s c ) = ^^F(s a ), (3.29) 
defining the function of one variable 

F ( s ») = «TFT? — V / Sd+(S ° ds <* ( 2s <i + s -- 3s °) ( 3 - 3 °) 

8K J (S a ) Js d _( Sa ) 

■F(3s - s - Sd, S a , Sd). 

This will allow us to rewrite the integral equation for F(s a , Sf,, s c ) as an equivalent 
set of coupled integral equations for functions of one variable only. 

Note that the projections as given in eqs. ( |3.18| ), (|3.19| ) and ( |3.30[ ) are defined only 
for values of s a in the physical decay region 4 < s a < (m — l) 2 . Outside this range 
the argument of the square root in the Kacser function eq. ( |3.17| ) becomes negative 



and we have to give this funcion a meaning by appropriate analytic continuation. 
This was studied by Bronzan and Kacser |27], 28 1 whom we follow closely. Based 



on comparison with explicit expressions for the projection operation in perturbation 
theory and demanding that the general definition of the projection should reduce to 
these expressions, they find that the naive integrals have to be replaced by contour 
integrals in the complex s^-plane. The respective paths joining Sd-(s a ) with Sd+(s a ) 
must avoid the real axis for 4 < Sd < oo as the amplitude F has a cut there. Where 
necessary, Sd~(s a ) and Sd+(s a ) are taken infinitesimally above or below the real axis 
according to the prescription obtained by replacing m 2 —>■ m 2 + id, 5 — * 0+ for real 
s a - Hence, we redefine the projections as the contour integrals 

1 f s d+{Sa) 

F ( S a) = nrs/ \ 9 ds d F(s a , S d , 3s - S a - S d ) , 

IK {S a ) Js d _(s a ) 
1 f s d+{Sa) 

F ( s °-) = 7^71 — \ f ds d F(s d ,3s - s a - Sd, s a ), (3.31) 
2K{s a ) Js d _( Sa ) 

3 f s d+(Sa) 

F ( S «) = or/V \ f s d (2Sd + S a - 3S ) 

■F(3s - s a - s d , s a ,s d ) 

and introduce the function H 



H{s) 



\{s - A){s - (m - l) 2 )(s - (m + l) 2 ) 



(3.32) 



which is well defined for all s different from zero and allows us to express the nec- 
essary analytic continuation of the Kacser function K in a simple way. 

Next we turn to the description of the different paths in the Sd-plane belonging 
to the four cases we have to discuss if evaluating the aforementioned definition of 
the projection operation. We must distinguish four cases. 
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i) 4 < s a < m,2 ~ 1 : as s a is in the physical decay region K{s a ) coincides with H{s a ). 
(see Fig. la). The end points of the contour are lying infinitesimally above the real 
axis and are found to be 

s d+ (s a ) = ^^ + H(s a )+t5, (3.33) 

Sd~(s a ) = 3g ° 2 Sa - H(s a ) + iS. 

ii) m2 2 ~ 1 < s a < (m — l) 2 : s a is still in the physical decay region so that K(s a ) 
coincides with H(s a ) again. The point s a = m2 2 ~ 1 is important because when s a 
goes through it, the endpoint Sd-(s a ) continuously turns around the beginning of 
the cut at 4 and remains then infinitesimally below the real axis (see Fig. lb) such 
that only the other endpoint Sd+(s a ) of the contour remains infinitesimally above 
the real axis. The two end points are now found to be 

s d+ (s a ) = ^^ + H(s a )+tS, (3.34) 

s d -(s a ) = g ° 2 Sa - H(s a ) - i8. (3.35) 

iii) (m — l) 2 < s a < (m + l) 2 : as s a has left the physical decay region K(s a ) 
must be analytically continued and is defined to be K(s a ) = iH{s a ). The endpoint 
Sd-{s a ) turns now smoothly into the half plane below the real axis whereas the 
other endpoint Sa+{s a ) turns smoothly into the half plane above the real axis in a 
symmetric fashion such that 

Sd+( s a) = 3S ° 2 Sa +iH(s a ), s d „{s a ) = 3S ° 2 Sa -%H{s a ). (3.36) 

The most adequate contour for numerical purposes joining those two points follows 
just the path the end points move along if s a goes from (m — l) 2 to (m + l) 2 (see 
Fig. lc) and is thus parametrized by s itself as shown in eqns. ([3.36|) . 

iv) (m + l) 2 < s a : as s a is out of the physical decay region K(s a ) must again be 
analytically continued and is defined to be K(s a ) = —H(s a ) (see Fig. Id). The end 
points of the contour are lying in this case on the real axis 

s d+ (s a ) = 3g ° 2 Sa - H(s a ), s d ^(s a ) = 3g ° 2 Sa +H(s a ). (3.37) 

Note that at s a = (m — l) 2 the projections of a function with a cut have a 
singularity because K vanishes at this point whereas the integration contour still 
has finite length and the integral over the imaginary part of the function along 
the contour does not vanish. The actual form of this singularity for F(s a ) around 
s a = {m — l) 2 is then given by 

F(s a ) ~ ± -r " const. (3.38) 

{(m - l) 2 - s a )2 
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and for F(s a ) by 

F(s a ) ~ J- nr • const. (3.39) 

((m - l) 2 - s a )2 

where the constant is in fact the value of the contour integral for s a = (m — l) 2 . We 
will use this form later to analyse the consequences of this pole in the projection 
prescription. 

The projection operations occuring in the eqs. (|3.18| ), ( |3.19| ) and ( |3.30| ) are now 



given by well defined contour integrals and we use them in the further evaluation of 
the integral equation ( |3.1(J| ) for F. We thus insert the partial wave expansions of T 1 



in the expression ( |3.10|) for the scalar function F. Using the definitions (|3.18|) , (|3.19| ) 
and ( |3.30D for the different projections and including the subtraction polynomial as 
required by the discussion in section 4, we may recast F in the following way 

F(s a ,s b ,s c ) = P(s a ,s b ,s c )+U(s a ) + V(s b ) + V(s c ) (3.40) 

+ (S c - S a ) • W( Sb ) + (S b - S a ) • W(s c ) 

where we defined the functions of one variable U, V, W 

U(s a ) = nil(ga " gt) f =rn % ~, (3.41) 

■{f°o*(s' a )F(s f a ) + l(f° *(s' a ) - f *(s' a ))F(s> a )}, 



VM _ ny^^ y } (3 , 2 ) 



7T 



and 

(s a - si)(s a - s 2 ) f°° ds' a 



W(s a ) 



3tt Ja (s' a - si)(s' a - s 2 )(s' a - s a -ie) 

■fl\s' a )F{s' a ) (3.43) 



+ perm., 



where perm, denotes the symmetrization of W in s\, S2, S3. Note that, unlike U and 
V, the function W is subtracted only twice. As ChPT yields the amplitude F only 
to the accuracy p 4 and as W is multiplied by factors of 0(p 2 ) we can not make three 
subtractions here. On the other hand, convergence of the integral is still ensured as 
the p-wave projection brings an extra factor of \ coming from the Kacser function 
K in the integrand. 

Numerically it is easier to deal with three integral equations for three functions of 
one variable than with one such equation for a function of three variables. We thus 
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have to transform the content of the subtracted dispersion representation eq. p.lOj) 



for F by taking its s- and p-wave projections. Note that all projections of functions 
of one variable are either trivial or may be expressed as a "bar-projection" (see eq. 
(|3.180 , where each argument s c is understood to be expressed by the integration 



variable s b and constant terms by 

s c = 3s -s a - s b . (3.44) 

The projections become now 

F{s a ) = P(s a ) + U(s a ) + V b (s a ) + V c (s a ) 

F(s„ 



+ 


(S c ~ 


Sa)W b ( 




P(Sa) 


+ U b {s 


+ 


(S a ~ 


s b )W c ( 






3 




+ AK 2 



F(s a ) = P(s a ) + ———{(s b -s c )U c (s a ) + (s b -s c )V b (s a ) 



+ (s b - s c ) 2 ■ W(s a ) + (s b - s c )(s a - s c )W b (s a )} 

where the subscripts b, c denote the argument of the respective function under the 
projection integral. Note that certain terms vanish because s b — s c = 0. We next 
change arguments of the functions above to s b and omit in the following the sub- 
scripts making use of the relations 

A c (s a ) = A b (s a ) 



s b A c (s a ) = s c A b (s a ) (3.45) 
s c A c (s a ) = s b A b (s a ) 

where A denotes any of the functions U, V, W. Eqns. (|3.45| ) are thus rewritten as 

F(s a ) = P(s a ) + U(s a ) + 2V(s a ) 

+ 2{sc-s a )W{s a ) ) 
F{s a ) = P( Sa ) + U( s a ) + V(s a ) + V(s a ) 

+ Oa - S c )W(s a ), 

3 



F(s a ) = P(s a ) + AR2 , {(s c - s b )U (s a ) + (s b - s c )V (s a ) 



1 a) 



+ (S b - S c f ■ W( Sa ) + (S b - S c )(.S a - S C )W ( So )}. 

Finally, we insert eqs. ( [3.461 ) for the projections expressed in terms of U, V, W 
into the defining eq. ( ft.41J) for U, eq. ( |3.42| for V and (|3.43|) for W. This yields the 
coupled system 

tt( \ n?=i( s a — Si) f°° ds 
U{s a ) = — 



7T 



4 Ilj =1 (s' a - sj)^- s a -ie) 
14 



[TV.) + \P{B' a ) + U( S ' a ) + p( S > a ) 

+\v{s' a ) + \v{s' a ) + \{s c -s a )W (s' a )} (3.46) 

-/o 2 Va) [f^O + \Wa) + fv(0 

+^(0 + 1^-^^)]}, 



nLi(sa-Si 



7T 



n 



j=l( S a 



/<f(<) + g(Q + V(s' a ) 

+V(s f a ) + (s a -s c )W(s' a )} 



(3.47) 



(fa ~ glKga ~ gg) f 00 ^ 

3tt A (s' a - si) « - s 2 ) « - s - ie) 

•A 1 * (4) [P(s'a) + ^^{(sc - s 6 )^«) + (s b - s c )V(s' a ) (3.48) 

+ (s b -s c ) 2 ■ W(s' a ) + ( Sb -s c )(s a -s c )W «)}] 
+ perm. 

These are our final expressions which generalize the Khuri-Treiman equations. As 
we mentioned above there are in fact singularities in the projected functions occuring 
under the integrals. The question of whether this leads to problems in the recon- 
struction of the amplitude by the dispersion integrals will be analysed in section 
5. 



and 

W(s a ) = 



4 Subtractions 



The generalized Khuri-Treiman equations in (|3.40| - |3~43|) have been formulated with 



three subtractions, because we will use them in this form. However, as usual in dis- 
persion relation techniques, the number of subtractions is to a large extent arbitrary. 
The minimal number of subtractions is determined by the high energy behaviour of 
the amplitude. If this is known, one is of course still free to take more subtractions 
than necessary for convergence. With fewer subtractions, the dispersion integral 
is more weakly convergent and depends more strongly on the poorly known high 
energy behaviour of the scattering phase. Conversely, if many subtractions are used 
the integral is better behaved; however then we need a larger number of subtraction 
constants which in general are less well known. Thus, one must keep a delicate 
balance in order to optimally use the available data and their uncertainties. 
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The asymptotic behaviour of the rj — > 3ti amplitude may be indicated by Regge 
phenomenology. Since there is no pomeron exchange in t]tt — > Tin, the leading 
Regge trajectory is the one associated with the p; we therefore expect an asymptotic 
behaviour A ~ s l J 2 (s b = 0) for the 1=1 component of the amplitude in the t- 
channel. The 7 = 0,2 components are not affected by p-exchange and are expected 
to approach a constant asymptotically In order to ensure convergence, it is thus 
sufficient to subtract the dispersion relation once. 

However, our calculational scheme of dealing with nn final state interactions relies 
on elastic unitarity. In the 1 = s-wave channel the phase shift exhibits large 
inelasticities above 1 GeV. Therefore it is desirable that this region is unimportant 
in the evaluation of the dispersion integrals. It turns out that this requirement 
implies at least two subtractions. We will discuss the influence of the phase shifts 
above 1 GeV on our results in sect. 6. 

Besides the final state interactions we are attempting to control, there are also 
mass corrections to the rj 3n amplitude. These are also known at the one-loop 
level where the total amplitude is M 

s b s)- g ' 2 t^~<) x 
{T{s a ) + U(s a , s b , s c ) + V(s a ) + W{s a , s b , s c )} , (4.1) 

where according to the discussion in section 1, the normalization of the amplitude 
is given by the first numerator rather than by 

m\ = m 2 K o - m 2 K+ - m 2 + m 2 + (4.2) 

which would follow if Dashen's theorem is used to determine the electromagnetic 
mass difference of the kaons. The leading order expression T(s a ) is given in eq. 
(|1.1|) and the other term will be discussed below. U(s a , s b , s c ) denotes the unitarity 
corrections and V(s a ), W(s a , s b , s c ) are polynomials in s a ,s b ,s c resulting from tad- 
pole and tree graphs at order p 4 in the low energy expansion. U(s a , s b , s c ) contains 
the two-pion I = 0, 1 and 2 final state interactions. These rescattering graphs are 
primarily responsible for a large enhancement of the amplitude at next-to-leading 
order. However, as discussed in great detail in Ref. 0, only the sum of all terms 
in eq. ( |4.1p is a meaningful quantity: the magnitude of the final state interactions 
of the pions depends on the manner in which the unitarity corrections are split off 
from the rest. 

In this work we use the ChPT one-loop amplitude to obtain an improved estimate 
for the subtraction polynomial entering the generalized Khuri-Treiman equations. 
There are three important points here: i) As mentioned above the one-loop am- 
plitude contains already some final state interactions. In order to avoid double 
counting, these final state interactions have to be incorporated in a well defined 
manner, ii) The choice of the subtraction points: In principle, a dispersion relation 
can be subtracted at any point. However, since the amplitude is approximated by 
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a polynomial near these points, they should lie as far as possible from kinematical 
singularities. Furthermore, we expect that the low energy expansion is more reli- 
able, if the invariant mass of the interacting pion pair is as small as possible at the 
subtraction points. However, as soon as I = 2 (and 1=1) final state interactions 
are included, the problem depends on two variables s a and s&, i.e. the pions interact 
also in the Sb- and s c -channels. This is reflected by the fact that there are three 
dispersion integrals. They should all be subtracted at small momenta, but there is 
no point in the Mandelstam plane where all three variables s a , s& and s c are small, 
iii) The subtraction constants as determined from the existing one- loop calculation 
will have uncertainties due to yet unknown higher order effects. These uncertainties 
can be estimated by comparing the one-loop amplitude with the leading order term. 
A precise description of our estimate will be given in section 6. 

The solution to the first two problems is straightforward and motivated by the 
following observation. In chiral perturbation theory, the imaginary part of the r\ — > 
3tt one-loop amplitude is generated by the graphs where tttt, ttt], KK or rjr] rescatter. 
They are easily obtained from two-body unitarity provided the lowest order vertices 
are inserted in the unitarity relation. Given the absorptive parts, the real part at 
order p 4 can be reconstructed from a dispersion representation up to a second order 
polynomial in s a , Sb, s c . Due to the work of Gasser and Leutwyler 0, this polynomial 
is known in terms of m n , tuk, m v , F n , Fk and the low-energy constant L3. 

The Khuri-Treiman equations are of course more general than this dispersion rep- 
resentation of the 77 — > 37r amplitude to order p 4 . However, they can be matched 
onto the ChPT one-loop amplitude in the following sense. Consider an itera- 
tive solution of eqs. ( |3.40| - |3.43| ). In particular, use the current algebra result, 
F(s a , Sb, s c )^ = T(s a ), as a first approximation to be inserted on the right hand 
side of eq. ( |3.41| ) and the partial waves of the nn scattering amplitude as given by 
their low energy expansion, i.e. 

flW-eWtmSMnS^is) (4.3) 

where <5/' ChPT (s) are the expressions for the phaseshifts at leading order ChPT given 
in Ref . . By construction, the first step of this iteration yields exactly the imagi- 
nary parts of the ChPT one-loop amplitude which are due to ttti intermediate states 
[]. Moreover, the dispersive part of the Khuri-Treiman equations contains also all 
7T7r-threshold effects present in the ChPT one-loop amplitude. Therefore, after the 
first iteration step as described above, we have reproduced all singularities close to 
the physical region. We note already at this point that three subtractions will be 
necessary in order to make this procedure well defined. 

In the low energy region, the remainder of the ChPT one-loop amplitude can be 
expanded to a very good accuracy into a polynomial of second order in s a , s&, s c . 

5 The absorptive parts due to iri], KK or 7777 intermediate states were not included in the 
calculation of the imaginary part. Because of their high thresholds, we expect a small effect and 
will neglect them. 
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This last step then determines the subtraction polynomial of the generalized Khuri- 
Treiman equations. 

Before turning to the details of our subtraction prescription, we give a more 
complete discussion of the structure of the one- loop amplitude in Chiral Perturbation 
Theory. 



4.1 Structure of the ChPT one-loop amplitude 

We work with the reduced amplitude A defined by 



,2 



A^(s a , s b , s c ) = ^=|^ -A(s a , s b , s c ) (4.4) 

where A 1_loop (s a , s b , s c ) is given in eq. (|4.1|) . T, V and W are polynomials to second 
order in s a ,s b} s c . The aim is to rewrite y4 1_loop (s a , suitable dispersion 

integral. The singularities of the amplitude A are contained in the unitarity correc- 
tions U(s a , s b , s c ). We may split them up further by writing 

U(s a ,s b ,s c ) =U* i n sp (sa,s b ,s c ) +U iem (s a ,s b ,s c ), (4.5) 

with 

U^(s a , s b , s c ) = ^A ^(s a ) [3T(s a ) + T(s b ) + T(s c )\ 

+iA 1)7r7r (s 6 ) [T{s a ) - T(s c )} + -Ai i7r7r (s c ) [T(s a ) - T(s b )} 
+2 A 2i7r7r (> 6 ) [T(s a ) + T(s c )) + ~A 2)7r7r (> c ) [T(s a ) + T{s b )) 
--A2,^(s a )[s b (t)+T(s c )} (4.6) 

and 

„2 -oo r/,ChPT/ A 

= 7 L" w-.-ly (4 ' 7) 

The phaseshifts S I,ChPT are the expressions to leading order in ChPT given in Ref. 
[§. contains the singularities due to tttt intermediate states. The remainder, 

[/ rem , is polynomial except for singularities due to nr], KK or r\r\ intermediate states. 
More explicitly, it is given by 

U^ m (s a ,s b ,s c ) = ^A (s a )[3T(s a )+T(t) + T(s c )} 

+^A 1 ( S6 ) [T(s a ) - T(s c )} + ^Ai(s c ) [T(s a ) - T(s b )} 

+ ^A 2 (s b ) [T(s a ) + T(s c )} + X - A 2 (s c ) [T(s a ) + T(s b )} 
+A 3 (s a ) (4.8) 
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with 



A (s) 
A 1 (s) 



1 

F 2 



(2fc* 



1927T 2 



^ {-2m^_ + (Aw + g^) s } 
(3s 



Am 2 K ) 



KK 



» + ^ j: 



3F 2 n 



A 3 (Sa) 



F 2 



— -_(T( S6 )+T( Sc )) 



96tt 2 



s a (3s a - Ami] 
AF 2 (m 2 — m 2 



"•^K"i<:( s a) + 



m 2 (3s a -Ami) 
3F 2 (m 2 -Mvr 2 ) 



ml 



3s a (3s a 



2F 2 rmK <*> gjp2( 



Jkk{so) 



TTT} 
1 



(4.9) 



The functions fcpp, J P q, Jpg and M PQ are defined in Ref. and are not displayed 
here explicitely; when not given, their arguments are as on the left hand sides of 
the equations. U^ p (s a , s b , s c ) will be reproduced by iterating the Khuri-Treiman 
equations as discussed before and shown in more detail below. In fact, it is uniquely 
determined by unitarity. The rest of the amplitude, 

A iem (s a ,s b ,s c ) = A(s a ,s b ,s c )-U^ p (s a ,s b ,s c ) 

= T(s a )+V(s a ) + W(s a ,s b ,s c )+U rcm (s a ,s b ,s c ) (4.10) 

is smooth in the physical region and it is expanded in a second order polynomial in 
s a , s b , s c around the center of the Dalitz plot, s a = s b = s c = s = \{m 2 + 3ml): 



Arem(s a , s b , s c ) = a + (3(s a - s ) + j(s a - s ) 2 + 5(s b - s c ) 



(4.11) 



This form is the most general second order polynomial if the constraint s a + s b + s c = 
3so = |(m 2 +3m 2 ) and bose symmetry are satisfied. The expansion of U rcm is tedious 
and its analytic form will not be displayed here. T(s a ),V(s a ) and W(s a , s b , s c ) are 
already polynomial and their contribution to a, 5 is 



a 



T(s )+V(s ) + W(s ,s Q ,s ) 
3 

{1 + a x + 3a 2 (m 2 - m 2 n ) + a 3 (9m 2 - ml] 



(m 2 - ml] 



4m 2 



[mi 



7 
5 



6 



F 2 (m 2 - ml] 



F 2 (m 2 



mi 



d2)}~ 



12s 



FZ(m 2 v 



mi 



mi 



(4.12) 
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The constants a%, 0,2, 03, d\ and d<i are tabulated in Ref. ||. We note here that 
the coefficients j3, 7 and 8 depend on the low energy constant L 3 which is phe- 
nomenologically not known very accurately Q. Its value extracted from decays 

is m 

L 3 = (-3.62 ± 1.31) ■ 10" 3 . (4.13) 
The error bar on L 3 leads to corresponding error bars on j3, 7 and 8: 

Ai. a = 

1,3 

12s 

^ = - q K ^) W - *>.76GeV- 



7T V 77 7T/ 



-1 



= - FZ im? ,-,n V) (AL ^* 1MGeV (414) 

We shall comment on the importance of this uncertainty in the determination of L 3 
in section 6 where we discuss the phenomenological implications of our results. 

The information needed to fix the subtraction constants of the generalized Khuri- 
Treiman equations is contained in coefficients a, . . ., 8, eq. ( |i.l2|) , as well as in the 
corresponding expressions coming from the expansion of U rem . 



4.2 Fixing the subtraction constants 

We consider the generalized Khuri-Treiman equations in the three times subtracted 
form given in eqs. ( |3.40| )-( |3T4"3|) Q The si, s 2) ... w 2 are finite subtraction points, 
fi = sin dfe^i is the partial wave for nir scattering with angular momentum I in the 
isospin channel / and the precise definition of F and F has been given in section 3. 
P is a polynomial of second order in s a , s&, s c with the same form as eq. ( |4.11|) : 

P{s a , s b , s c ) = a + (3s a + 'jsl + 5(s b - s c ) 2 . (4.15) 

The reason for using three subtractions will become clear soon. Note that the 
polynomial P is not equal to A Icm because the latter refers to a twice subtracted 
dispersion relation (see eq. (|4.7| )) while P corresponds to three subtractions. 

Now we iterate eq. ( j3 .41]) by inserting on the right hand side latexthe current 
algebra expression F(s a , s c ) = T(s a ). Moreover, we replace // by its low energy 
expansion 

//(a.)«^ ChPT (a-)- (4-16) 

We observe that this reproduces the absorptive parts of the ChPT one-loop ampli- 
tude due to mi intermediate states. For isospin I = 0, 2 this is trivial because T(s a ) 

6 Note that L3 is scale independent. 

7 Note the functions U, V, W are not identical to the functions U, etc. introduced above. 
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is linear in s a and hence T(sb) + T(s c ) is again a first order polynomial in s a . For 
the 1 = 1 p-wave nir intermediate state the absorptive part is in general given by 

lmF(s a , s b , s c ) lsl = e i6 ^ sm8\{s b ) ^^ F^{s h ) + (t <-> s c ) (4.17) 

where K(sb) is the Kacser function and F^^Sf,) is the p-wave projection in the s&- 
channel of the function F(s a , s&, s c ) defined in section 3. Setting T(s a ) = a + frs a , 
the p-wave projection in the s fc -channel is calculated to be 

2?" (s b ) = -\K{s b )b. (4.18) 

Thus the first iteration with the approximations given above yields 

lmF(s aj s b , s c )| = # ChPT ~(r(s a ) - T(s c )) + ( S6 <-> s c ). (4.19) 

This is indeed the result of the one-loop calculation obtained in Ref. ||. 

Turning now to the dispersive part, we shall show how the real part of W^ sp , eq. 
(|4.6|) , is reproduced in the first iteration of the generalized Khuri-Treiman equations. 
In order to make the argument more transparent, let us consider first the simpler 
case where all subtractions are taken at zero. The term U(s a ) on the right hand 
side of eq. ( |3.40| ) then reads, after the first iteration 



,3 f00 rO.ChPT, ,w~ La _ ,2,ChPT , r n 



with 



and 



5 i 2 i 

-a + bso, 3 

\a + bs , h= -\ h ( A2V > 



3 s 



ml — ml ' ml — ml 



(4.22) 

"Tj 7T ""1J 7T 

In eq. ( |4.2U|) , P denotes the principal value of the integral; the combinations in eq. 
(|4.21|) simply stand for the appropriate isospin projections. Using the decomposition 

~ P / ds ^TT^ ^ = ( s - Z ^~ P / ds ,2 ( , ■ \ + Z ~ P / s — ~rT~ 

7r J s 6 {s ■ — s — ie) n J s z (s •' — s — le) n J s a 

(4.23) 
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eq. ( f4.20| ) can be brought into the form 



ReU iter (s n 



T(s a ) + -(T(s b )+T(s c )) 



ReA 0i7r7r (s a ) 



0,ChPT, 
\ 



-,/3 



-a—P I ds 

IT J 

~[T(s a ) +T(s c )}ReA 2 ^(s a ) 



s 2 f 
+a — P / ds 
ix J 



A 



2,ChPT, 
I 



-./3 



(4.24) 



We thus have transformed Re£7 lter (s) into the first and last term of IA^ V plus two 
polynomial terms. In other words, eq. ( [4.23| ) serves to transform a three times 
subtracted relation into a twice subtracted one. Likewise, we may decompose 



Re(V(s b )+V(s c )) 



iter 



~[T(s a )+T(s c )]ReA 2)7r7r (s 6 ) 



with 



-a—P I ds 



a = a + 2^ s °- 



2,ChPT, 



,/3 



+ (s b <-> S c 



(4.25) 



(4.26) 



Finally, using ( J4.18|) and the definition of F(s'), eq. ( |3.38| ), we see that the 1 = 1 
contribution in the generalized Khuri-Treiman equations yields precisely the terms 
in proportion to Ai j7r7r in eq. ( [4.61) : 

Re ((s a - s c )W(s b ) + (s a - s b )W(s c )) itCT = X - [T(s a ) - T(s c )] ReA 1)7r7r (s 6 )+(s 6 ^ s c ). 

(4.27) 

We see that linking the dispersion technique to the one-loop ChPT calculation al- 
lows to subtract U,V and W separately avoiding the problem (mentioned previously) 
that these points may lie outside the physical region. 

We are now ready to determine the subtraction polynomial P of the triple sub- 
tracted Khuri-Treiman equations. After the first iteration, the amplitude can be 
written as 



F(s a , s b , s c ) itei = P(s a , s b , s c ) + U^ sp (s a , s b , s c ) + R(s a , s b , s c ). 
If all subtractions are taken at zero, R is given by eqs. (|4.24 )- ( |4.25[ ), i.e. 



R(Say "Sb; 



-a^P I d/° 

7T 



O.ChPT/ n 



-,/3 



\ 7T 7T 



P / ds 



A%> ChPT (s> 



-./3 



(4.28) 



(4.29) 
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Requiring the matching of the first iteration of the generalized Khuri-Treiman equa- 
tions with the one-loop amplitude from ChPT we obtain 



= a + (3(s a - s ) + 7(s a - s ) 2 + 5(s b - s c ) 2 - R(s a , s b , s c ) 

where the subscript "expand" means a Taylor expansion up to second order in 
s a , Sb, Sc', in fact the expression in the curly brackets is just A Tem — R. 

Eq. ( |4.30| ) and its generalization for finite subtractions points described below is 
the main result of this section. It describes our method of using the ChPT one-loop 
amplitude (see however below) to fix the subtraction constants of the generalized 
Khuri-Treiman equations. In order to make the procedure transparent, that is show 
explicitely how the higher corrections unitarize the amplitude, we were forced to 
subtract the Khuri-Treiman equations three times. Otherwise the integrals occuring 
in R(s a , Sb, s c ) would not converge. By construction, the first step in an iteration of 
the Khuri-Treiman equations reproduces the ChPT one-loop result in the physical 
decay region to very good accuracy. The numerical solution of the generalized 
Khuri-Treiman equations as attempted in section 5 can therefore be interpreted as 
a correction on top of the ChPT one-loop amplitude due to all possible two-body 
final state interactions of pions. 

There remains the important question of higher order corrections to the one-loop 
results which may shift the subtraction constants substantially. We will discuss 
adress this issue in section 6. Here we only note that also in this case the above 
framework of three subtractions can be employed if the necessary modifications are 
made. 

4.3 Finite subtraction points 

In the preceeding subsection all three subtractions have been taken at zero for the 
sake of simplicity. Here we shall describe the modifications which arise if the sub- 
tractions are taken at finite values. We do this for two reasons: Since the dispersion 
relations do not fix the subtractions points, varying them over a certain domain 
gives an estimate of the error of our procedure (see also the discussion in section 6). 
Second, it is numerically favorable to have subtraction points which do not coincide 
when we solve the Khuri-Treiman equations iteratively. 

Introducing now finite subtraction points and performing essentially the same 
steps as before, we find that only the function R(s a , Sb, s c ) is modified. It explicitly 
depends on subtraction points Sj,t>j, i=l,2,3 and w i} i=l,2 and is given as 





expand 



(4.30) 




+bR i2) (s b ; vi, v 2 , v 3 ; - ~ ) + W? (2) (s c ; Vi,v 2 , v 3 ; - y ) 



b , , UJ b 




(4.31) 
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The explicit form of functions I = 0, 1, 2 is given in Appendix D. 

The subtraction polynomial is then again given by eq. ( f4.30|) with R taken from 
(|4.31|) . It depends explicitly on the subtraction points Sj, t>j and Wi. However, 



by construction this dependence on the subtraction points is counterbalanced by a 
corresponding dependence of the dispersive part of the generalized Khuri-Treiman 
equations, provided only the first step of the iteration described above is performed. 
The procedure of fixing subtraction constants is in this sense independent on the 
choice of the subtraction points. Up to the first iteration, the ChPT one-loop ampli- 
tude is reproduced, in the physical region, for any value si, w 2 in the low energy 
regime. Beyond the first iteration this is no longer true. 

Finally we give numerical values for coefficients a, 5 of the subtrac- 
tion polynomial P. We consider the sets of subtraction points SP = 
Si = v\ = Wi, s 2 = v 2 = w 2 , s 3 = t> 3 displayed in Table 1. The input parameters 
we use are m n = 140MeV, m v = 549MeV, = 92.4MeV, F K = 114MeV, 
L 3 = (—3.62 ± 1.31)lCr 3 . The calculated coefficients a, S for these sets of 
subtraction points are given in Table 2. 

The dependence on the choice of subtraction points is substantial and we shall 
discuss how it propagates into the final numerical solution in section 6. Note that 
error bars due to the uncertainty in the determination of L 3 are correlated. For 
instance, the value of P(s a , Sb, s c ) at the center of the Dalitz plot does not depend 
on L 3 - the errors for a, 6 cancel for this quantity. Also it appears that for 
subtraction points SP ~ 0, 7 is rather small. In the physical region, the term in 
proportion to 7 contributes only « 2% to the subtraction polynomial. The error bar 
on 7 due to the uncertainty in L 3 is therefore phenomenologically less important. 



4.4 Comments 

We have swept over several subtle points rather briskly and would like to come back 
to them. 



1. The polynomial P depends on four constants. On the other hand, we have 
subtracted the three functions U, V and W separately; this implies 3+3+2 = 8 
constants. It is clear from the form of P that there are only four physically 
relevant parameters and that there are redundant parameters. The reason 
for this is found in a general invariance of the dispersion relations: they fix a 
function only up to a polynomial; therefore, we can 'shift' certain constants 
from one function to the other. For instance, we can redefine in this way two 
constants in W and two in V (the constant and the linear terms). Thus, only 
four constants remain. 

2. The asymptotic behaviour of the amplitude has been discussed only briefly. 
Taking three subtractions we want to make sure that the actual behaviour in 
the asymptotic region does not matter for the solution in the physical decay 
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region. There is however an other subtlety which we should mention. Math- 
ematically speaking, the solution to the integral equations ( |3.40| - j3l4"5| ) is not 



completely determined by specifying the asymptotic behaviour of the ampli- 
tude. The number of parameters (extracted from ChPT) to be input depends 
also on the asymptotic behaviour of the 7T7r-phase shift. Q Let us consider the 
twice subtracted relation for defmiteness and omit V and W. The solution to 
this problem can be written in the form 

F(s) = P(s) + <P 10 (s) + Ms) (4-32) 

where <Pi and ip are defined in eqs. (|5~8| - |5.13|) . Now, if the phase goes to zero 
for large s, the Omnes factor is constant; correspondingly, the amplitude rises 
linearly with s. On the other hand, if the phase approaches tt, the Omnes factor 
decays like s -1 and the amplitude tends towards a constant. Stated differently, 
a given asymptotic behaviour of the amplitude, say ~ const., requires one or 
two subtractions for the two limiting cases of the phase shift respectively. This 
ambiguity in the phase is reflected in an additional term 

T7r\( s ~ - s 2){s - s 3 ), Q(s) polynomial, (4.33) 
D(s) 



which can be added to the approximate solution $ w , c.f. eq. ( |5.28|) . It intro- 



duces an uncertainty in the final solution which should be well controlled since 
the assumptions entering the analysis are not rigorous. We have estimated this 
uncertainty by modelling the phase in the high energy region. The details as 
well as the resulting error bars on our final answer are given in section 6. 



5 The iteration 

With reasonable assumptions we now build up an iterative scheme for the numerical 
solution of the set of coupled equations for the projected amplitudes U, V, W. It 
consists of two distinct kinds of iterations. The first one accumulates contributions 
to V and W for a corresponding part of U while the second yields the contributions 
to U itself. 

The first type of iteration relies on the assumption that to lowest order only U 
contributes, but not V and W. It is motivated by the fact that pion-pion scattering 
at low enegies is dominated by the s-wave, 1 = channel as the s-wave, 1 = 2 
contribution is suppressed by a much smaller phase shift f$ -C /q- The 1 = 1 
contribution comes from a p-wave scattering and therefore is also suppressed. We 
may thus expect that in comparison to V and W, U will still yield the dominant 
contribution to F in our refined analysis. 

8 We are indebted to H. Leutw yle r for pointing out this to us. A more complete discussion of 
these issues can be found in Ref. 1122(1 ■ 
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To organize this first iteration we take the terms with the phase shifts /q , f\ as 
perturbations and introduce accordingly a counting parameter A. /q is taken to be 
of 0(A°) whereas 

fS - A/ 2 (5.1) 
fi - ^fi 

are treated as small parameters. In principle, one could introduce different counting 
parameters for the two isospin channels, but for sufficently many iterations in A, 
this is irrelevant. Expansion of U, V and W in a series in A yields 

oo 

U{s a ) = £A fc f4(s a ) 



oo 

V{s a ) = £A fe \4(s a ) (5.2) 
o 

oo 

W{ Sa ) = X> fc W*(s a ) 



where we set by assumption Vo(s a ) = Wo(s a ) = 0. We insert these expansions in 
the integral equations eqs. (|3.46|) , ( |3.47|) and (|3.48| ) and obtain by equating equal 
powers of A 



U k (s a ) 



Yli=l( s a-Si) r ds' a 



71 



4 nj=M- s j)( s *- s a-i£) 



■ifSVa) {{P(s'a) + lP(s' a )}5 k0 + U k ( S ' a ) + 2 -U k { S ' a ) 

+lv k (s' a ) + ^V k (s' a ) + ^ Sc -s a )W k (s' a )} (5.3) 
2^ 



V,(s n 



+ -V k _ 1 (s' a ) + -(s a -s c )W k ^(s' a )}}, 



U! =1 {s a ~ Si) f°° ds' a 



and 



W k (s a ) 



k h Uj=M - s i)( s 'a ~s a - ie) 
•/o*(4) [P(s' a )5 kl + U k - 1 (s' a ) + V k ^(s' a ) (5.4) 
+V k . 1 (s' a ) + (s a -s c )W k . 1 (s' a )] 

{s a - sx)(s a - s 2 ) f°° ds' a 



3tt Ji (s' a - si)«- s 2 )(s' a - s a -ie) 

■flVa) [P(s'a)hl (5.5) 
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+ 4X2( y ) {0c - S b )U k -! (s'J + (S b - S c )Vfc_i (VJ 

+ - s c y ■ W k - X (s' a ) + (s b -s c )(s a -s c )W k ^ (s' a )}} 
+perm. 

The structure of the iterations becomes clear if we rewrite the equation for U k in 
the form 

TtLi(sa-Si) r°° ds' a 



U k (s a ) = <P ok (s a ) + [U=l[Sa Sz) I 

7T J 4 



4 U 3 j=i(s' a - Sj)(s' a - s a -ie) 
•ffWa) PM + \Uk{s' a )\ (5.6) 



introducing the combination 

_ Yll = x{Sa-Si) [°° ds 'a 



k J± Uj=i(s' a - Sj)(s' a - s a -ie) 
4- 



-/o*«) |P(0«M + \Uk-i{s' a ) + \v k ^{s' a ) 
2^ 



+-V k ^(s' a ) + -(s a -s c )W k „ 1 (s' a )}}. 

For given C4_i, I4_i, W k -\ and their respective projections we may use eqs. (|5^ )- 
( |5.5| ) to determine V k and W k . This allows to compute the subtraction function 
<P ok as given in eq. ( |5.7|) . If we can also solve eq. (|5.3|) for U k , the iteration step 
is complete and we may perform it once again for k instead of k — 1. Note that 
the lowest contributions Uo, Vq and Wo are completely determined by the different 
projections of the subtraction polynomial P containing the input information of 
ChPT on the decay amplitude. 

As the process converges rapidly, we will have to perform a small number of 
iterations, typically four or five, to obtain a precision below one part in thousand. 
Of course we take now A = 1 as it was introduced as a counting parameter only and 
may sum the different contributions to obtain the fully iterated U, V and W. 

We turn to the second kind of iteration which yields the yet undetermined U k 



from eq. ( |5.3| ). We thereby follow the discussion of Bronzan |29| and Neveu and 
Scherk []3l| extending it to our case. As we are dealing with one variable only we 
omit the usual subscripts here. 

If in the equation (|5\6| ) for U k the second term in the square bracket were absent, 
we would have to deal with an Omnes type of integral equation describing ordinary 
two pion-pion rescattering without rediffusion terms coming from the third pion in 
the final state. In our case making use of solvable Omnes equations we are able 
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to recast the equation for U k in a form which is now accessible for a second kind 
of iteration 31 1. To this end we introduce two auxiliary functions, the first of 
which is defined by 



ds' 



4 m =1 {s' ~ Sj ){ S ' - S - IS) 



(5i 



This equation is indeed solvable in terms of <Pok |36| and yields 



^ofc(s) 



ttD( S+ ) 7 4 n?= 1 (s'-s j )(s , -s-ie) 



(5.9) 



where s + = s + is and D(s) is just the Omnes function corresponding to the phase 

(5.10) 



shift <5° 



D{z) = e * J * «'-* 
The second auxilary function ^ is then defined by 



2— 

U h (s) + -U k (s) = $ lk {s) + Vfc(s). 

Using this definition twice and in eq. (|5.8|) , we obtain an equation for ^ 

(is' 



(5.11) 



3 7T 



4 YlU( s '- s i)( s '- s - i£ ) 



■f^(s')Ms) (5.12) 
which is again of the Omnes type and has a solution in terms of the unknown U k 

ds' 



fo(s')D(s' + )-U k (s'). 



Eqn. (|5.11|) is next rewritten as 



U k (s) = <Pi k (s)+M s ) - ^U k {s). 



(5.13) 



(5.14) 



But the difference of the last two terms is just given by the integral on the r.h.s. of 
eq. ( |5.13| ) so that we finally obtain the desired form of the equation for U k [JT 



UJs) 



ttD(s + ) Ja Yl* =1 (s'- Sj )(s>-s-i£) 
fo(s')D(s' + ) 2 -U k (s'). 



(5.15) 
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@ik is already determined in eq. ( |5.9| ) as an Omnes inversion. We may now account 
for the second term in eq. ( 5.15Q by iteration, noting that at every step one multiplies 
with the phase shift /q. Introducing the counting parameter \i 

/o°^M/o° (5-16) 

we expand U k in a series in /i 

oo 

U k (s)=Y,V m - 1 $ mk (s) (5.17) 
1 

and obtain for m > 1 

■f°(s')D(s' + )p {m _ 1)k (s'). (5.18) 

Again a small number of iterations is sufficient as the process here converges 
rapidly, too. Taking now \i = 1, we sum up the different contributions to obtain the 
full iterated U k . 

Next, we recast the main contribution $ik to Uk in a more convenient form. With 
its use we may give $io explicitly in terms of some projections of the subtraction 
polynomial P and the Omnes function. 

We start with the dispersion representation for $ fc as given in eq. ( |5.7| ) 



= ^ r ~, I - Ms") (5-19) 



7T 



4 nU( s "- s i)( s "- s - i£ ) 



where the function A is given by 



A(s) = f° *(s)[{P(s) + 2 -P(s)}6 k0 + 2 -Vk(s) 

4- 



+ ^V k (s) + -(s c -s a )Wk(s)} (5.20) 
-f^(s) [|p(s)(J fcl + ~Uk-i{s) + ^fc-i(s) 



+-y fc _i(s) + -(s a - s c )W k -i (s)} 

and insert it into the dispersion integral for as displayed in eq. ([5.9|). First, we 
focus only on the second term which becomes 

UUis-s,) r A{s") 



nD(s + ) h nj=i( s "- s i) 

±r ds > fS(*'M*' + ) i52V 

7T h {s" - s' - ie){s' - s - ie)' y ' ' 
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As the discontinuity of the Omnes function across its cut is given by 
D(s + ) - D(s_) = -2iD(s + )f%(s) = -2iD(s_)f°*(s) 
we may represent it as 



D{z) 



1 r 00 

7T J 4 



MOW 

s' — z 



ds 



(5.22) 



(5.23) 



With this representation, the s'-integration in eq. (|5.21|) may be performed and 
yields the difference 

■\D{s"_) - D{s + )} 



S" — S — IE 

so that the expression ( 5.21]) becomes now 



m=i(s-sj 



ds" 



D{s"_)A{s") 



(5.24) 



(5.25) 



Inserting this result in eq. (|5.19|) $ofc cancels and one obtains 



^ifc(s) 



ds' 



7TD(S + ) h U.J=l( S ' - S j)( S ' - S - i£ ) 



2 ~ 



D(s'_){fS*(s') [{P(s') + -P(s')}5 k0 + -V k {s') 



+-V k (s') + -(s c -s a )W k (s')\ 
-f**(s>) £p( S ')6 kl + 2 -U k ^{s') + lv k ^{s') 



(5.26) 



+ -W 1 ( S / ) + ~(s a ~ 8 e )W k -! { S ')}}. 

Up to now we only relied on the dispersion representation of #ofe such that the 
discussion is valid for general k. From now on it is assumed that D(z) is the only 
function with a cut under the integral ( |5.26|) . This is correct only for k = 0. Using eq. 
( |5.22| ) to rewrite /q*(s)D(s_) we may recast the dispersion integral along the cut L 
as a contour integral along C starting at oo and going down to 4 lying infinitesimally 
below the cut L, turning there and going back to oo infinitesimally above the real 
axis such that 



$io (s) 



2inD(s, 



dz' 



c ir j= i(z' - Sj)(z' - s - is) 
2 ~ 



■D(z') [P(z') + -P{z')) 



(5.27) 



i=l 



Ds) 
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Q(s) is a polynomial which is not restricted as long as no asymptotic boundary 
conditions are imposed. The form (|5.27 ) is easily evaluated with the help of the 
residue calculus and yields after a little algebra the final result 




l-D(s) 
1 D(s) 
- s 2 )(a - s 3 ) 



D(si)-l 
D(s) 

D(s 2 ) - 1 
D(s) 

D(s 3 ) ~ 1 
D(s) 



(si - s 2 )(si - S3) 
(s - si)(s - S3) 



(5.28) 



(s 2 - si)(s 2 - s 3 ) 
(s - si)(a - s 2 ) 



(s 3 - si)(s 3 - s 2 ) 



3 



+n( s -^) 



D(s)- 



i=l 



Note that it is of course possible to set two or more of the subtraction points s« 
in $10 equal. As one has then to deal with multiple poles the expression for the 
corresponding residue becomes more involved. 

We have built up an iterative scheme for the numerical solution of the set of 
coupled integral equations for the projected amplitudes U, V and W and now turn 
to a short description of the numerical aspects of the iterations. 

To obtain the decay amplitude with the inclusion of the final state interactions 
as described by the Khuri-Treiman equation we need to know the three functions 
U(z),V(z),W(z) only for physical values of z lying infinitesimally above the real 
axis, z = s + ie. But in both types of iterations described above, projections of U, V 
and W have to be computed. Since the corresponding projection integrals involve 
contours lying in the complex plane, the functions must be known in a certain 
domain around the real axis. 

As discussed in section 3 there is one s a -region, namely (m — l) 2 < s a < (m + l) 2 , 
in which the projection integral has to be performed along a contour joining two 
points lying in the upper and lower complex half planes respectively which must 
not intersect the real axis from 4 up to 00 as the integrands have a cut there. As 
they are analytic elsewhere we are free to choose those paths which minimize the 
number of necessary lattice points in the complex plane. The two possible choices 
join the respective end points as given in eq. ( |3.36|) along the path the end points 
themselves describe as s a increases from (m — l) 2 to (m + l) 2 . In a second s a - 
region with m 2 ~ 1 < s a < (m — l) 2 we have to perform integrals along paths lying 
infinitesimally above and below the real axis such that we have to know all the 
functions involved for z = s — ie if 4<s<m + l. 

As we are doing numerical computations, we have to replace the connected do- 
mains just discussed by a lattice. The spacing of its points is adapted to the accu- 
racy required and to the possible occurence of numerically problematical points. Of 
course we will cut off the different integrals at some large value of the integration 
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variable which is determined on one hand by the phenomenological knowledge of 
the pion scattering phases involved and on the other by demanding that numerical 
results do not change sensibly if one varies the cut-off. The corresponding cut-off 
for negative values of s a is then automatically fixed by eqs. (|3.36|) . As a result we 
have to know all the functions on a lattice of points lying in the complex plane as 
displayed in Fig. 2. The actual numbers for the different cut-offs are given in section 
6 and their impact on our results is discussed there. 



6 Numerical Results 

We consider the decay function in the form (see eq. (|3.40| )) 



F(s a , s b , s c ) = P(s a , s b , s c ) + U (s a ) + V(s b ) + V(s c ) 

+ (s a - s c )W(s b ) + (s a - s b )W(s c ). (6.1) 

Here, s a , s b and s c are the kinematical variables defined earlier. P is a polynomial 
and U, V and W correspond to the I = 0, 2 and 1 pion rescattering channels. They 
satisfy dispersion relations, the generalized Khuri-Treiman equations, which describe 
the unitary corrections to the 77 — > 37r decays (see eqs. ( p.46|J3.47| , |3.48|) ). We solve 



the dispersion relation iteratively, starting from an aproximative solution obtained 
by Neveu and Scherk. 

Using the existing one-loop results of Gasser and Leutwyler as basis and casting 
their unitary corrections into a dispersion relation, we are lead to a three times 
subtracted dispersion relation for the three amplitudes U, V, W. However, in order 
to assess more clearly the relevance of our subtraction procedure, it is useful to 
consider also the twice subtracted case where the subtraction polynomial is simply 
taken as the current algebra result for the amplitude. We will therefore discuss both 
forms of the dispersion relations, starting with the three times subtracted one. 



6.1 Uncertainties 

We can distinguish several sources of uncertainties. The first is associated with 
the ambiguities inherent in the dispersion relations. These relations do not fix 
the subtraction points, except of course for the requirement that they should lie 
in a region where the approximate theory makes sense. In addition, there is the 
ambiguity to the solution of the Khuri-Treiman equations related to the asymptotic 
behaviour of the 7T7r-phase shift as discussed in sect. 4. Then, there are the errors 
in the input parameters, that is in the values of the subtraction constants. These 
uncertainties are of course not completely disconnected; for instance the one-loop 
result from which we determine the subtraction constants may be rather accurate 
for some values of the subtraction points but not for others, and so a wrong choice 
of the subtraction point may underestimate the error completely. And finally, there 
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are technical problems, such as the convergence of the iteration and the numerical 
integration in the complex plane. 

1. Subtraction points 

We consider first the variation in the subtraction points. In section 4, we con- 
structed the second order subtraction polynomial which was obtained by writing 
the one-loop result in the form of a three times subtracted dispersion relation. This 
polynomial then serves as the starting point for the complete dispersion relation. As 
discussed in section 4, the choice of the three subtraction points sx, «2, S3 for U and 
those for V and W, which will be denoted collectively by SP, affects the subtrac- 
tion polynomial considerably and consequently the complete amplitude. Whereas 
in chiral perturbation theory and in view of the treatment in section 4 it seems most 
natural to subtract at small SP (at s, ~ 0), varying the points gives us a feeling for 
the error. We have therefore calculated some important quantities for several SP; 
the results are shown in Tables 3 and 4. As subtraction points we have taken the 
sets in Table 1. 

In Table 3 we give the splitting of the one- loop amplitude (0) in the center of 
the Dalitz plot into a polynomial and a dispersive part, as described in eq. Q4.30|) . 
Furthermore, we list the first approximation, $10, and our final result. Whereas, by 
construction of the subtraction polynomial, the one-loop amplitude at the center of 
the Dalitz plot is reproduced for all SP, the relative contributions of the two pieces 
vary strongly. However, as long as the subtraction points are not too close to the 
two-pion threshold, the total amplitude given in the last column of Table 3 is rather 
stable. We therefore consider it most natural to select a small subtraction point. 
Excluding SP = 3, we obtain A tot (s , s , s ) = 1.57 ± 0.12 + «(0.41 ± 0.03). 

Table 4 contains the physical observables for the various subtraction points. The 
rates in the second last two columns are calculated from eq. (|1.9|) , e.g. with the 
normalization of the amplitude given by ml- We will discuss the relevance to the 
quark masses in the next section. 

As may be expected, large subtraction points (near the physical region) yield 
smaller rates because there the full result is nearer to the one-loop result which is 
too small. On the other hand, large negative SP do not arbitrarily increase the 
rates; rather, these remain remarkably stable. Thus it is not possible to (artificially) 
enhance the rates by playing with the subtraction points. The uncertainty of the 
subtraction point can be substantially reduced, if the Dalitz plot slopes a, b, c (to be 
discussed below) could be measured with two percent accuracy; This would allow 
to narrow down the error bar on T to 10 % or 16 eV. At present, the experimental 
values (which admittedly have larger errors) prefer small SP. 

The values of the physical quantities which correspond to SP = will be taken 
as the central values. 

2. Errors in the subtraction constants 
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Without a next order calculation, the error on the subtraction constants in the 
polynomial, in particular on 7 and 5 which vanish at tree level cannot be truly 
assessed. A reasonable estimate of the errors can be obtained by the following 
observation by Anisovich and Leutwyler [22]. The current algebra result for the 
?7 — ^ 37r amplitude has an Adler zero (a value of the kinematic variables where the 
amplitude vanish), which for finite quark masses is shifted to s = sa = \ rr ^- This 
can be easily seen from the tree level result; obviously, the Adler zeroes lie on a 
straight line. The one-loop amplitude shares this feature (with a slightly shifted 
sa), if i n addition we fix s a = s c = sa (or s a = Sj = sa)- Moreover, along the line 
s a = s c the slope 4r(s a = s c = sa) is practically unchanged compared to the leading 



order expression 



This suggests that the amplitude along the line s a = s c 



near sa is very stable against corrections and that this kinematical point is therefore 
well suited as subtraction point. Note that it is crucial to stay on the line s a = s c . 
Consider for example the subtraction polynomial in the form corresponding to an 
expansion around the point s a — 0, s& = s c , 



P = a + (3s a + 7^a + S(s b - s c 



(6.2) 



The coefficients a, (3 in ( |6.2|) are very different from the current algebra expressions 



a 



CA 



-4my(m^ — m*), (3 = Zj(mi —m%). The amplitude is not stable against 
corrections along the line s& = s c and it is difficult to assign an error bar to a, ...5 
in ( |6.2|) as noted before. 

We therefore expand around s a 



s c = s A , i.e. 

, 7 , -2 



P = bz + cz + dw + ew + fwz, 



(6.3) 



where 



3s - 2s A - s b , 



w 



Due to Bose symmetry only three of the five constants in 
instance 



(6.4) 

are independent; for 



Numerically we find 



b-d 



6(s - sa) 

b-d 
2(s - s A ) 



+ c 



2c 



(6.5) 



c 
d 



4.87GeV" 2 

13.65GeV" 4 

11.88GeV -2 . 



(6.6) 



Here, only those contributions from the subtraction polynomial which come from 
A iem (see section 4.2) have been included. The remainder R(s a , s&, s c ) serves only to 
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bring the one-loop amplitude into the form of a dispersion relation and is cancelled 
once the first iteration step has been performed. As discussed above, the constants 
sa and b are known with very good accuracy. However, c receives contributions 
only at one-loop order and d is very sensitive to corrections. For the error estimate, 
we therefore assume that sa and b are given exactly by ( |6.6| ) whereas c and d are 
assigned a relative error of 25 %, typical of higher loop-corrections: 

Sc = 3.41GeV~ 4 

8d = 2.97GeV~ 2 . (6.7) 

We can now translate these errors into uncertainties of the constants a, S in the 
expansion ( |6.2| ) by comparing the two forms for P (again, the contributions of the 
remainder R(s a , s b , s c ) are not written): 

a = -1.02 ± 0.01 =F 0.34 

= (20.35 =F 0.72 ± 5.69)GeV~ 2 

7 = (-1.38 ± 13.65 T23.72)GeV" 4 

5 = (6.22 ± 0.0 ±2.63)GeV~ 2 . (6.8) 

The first and second error is induced by the uncertainty in c and d respectively. The 
errors induced by d are particularly large, however they are correlated and cancel 
completely for the polynomial (|6.2j ) evaluated at the center of the Dalitz plot. 

The induced error bar on the amplitude at s a = s c = Sq can be estimated by 
calculating the change in P(sq, sq, Sq). Only Sc matters of course and we find 

5A(s , s , s ) « 5P(s , s , s ) = 5c(s a - s A ) 2 = 0.12. (6.9) 

This is 8 % of the full amplitude corresponding to a 16 % error on the rate, i.e. 26 
eV if Q = QDashen is employed. The actual change of the full amplitude is somewhat 
larger and summarized in Table 7. 

The rate is thus very sensitive to the value of the constant c. If its higher order 
contributions turned out to be large and positive, using its one-loop value would 
seriously underestimate the rate. We therefore attempt a more quantitative estimate 
of the higher order corrections to c by considering a twice subtracted dispersion 



relation, following again Anisovich and Leutwyler p2| . We start by decomposing 
the amplitude into its isospin components, i.e. 

F(s a , s b , 8 e ) = M (s a ) + (s a - s c )Mi(s 6 ) + (s a - s b )Mi(s c ) 

+M 2 (s b ) + M 2 (s c ) - ^M 2 (s a ). (6.10) 

This decomposition is ambiguous for unitarity determines only the singular part of 
the amplitudes M , M\ and M 2 . A polynomial in s a can always be shifted between 
the isospin amplitudes, a feature which can be used to make the I = 1, 2 parts in F 
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small, throughout the physical region. Then, a dispersive analysis is performed for 
the function M only, which must satisfy the twice subtracted relation 



M (s a ) = -1- Lo + Po(s a - s A ) + {Sa SA? { £ 

D(s n ) IT J Is' 



Ur 

s A y(s' - s a -ie) 



s a -s A ) 2 f ds'D(s')f$(s%M (s>) \ 

(6.11) 

For simplicity, we have neglected the small pieces proportional to Mi, M2 under the 
integral. The two subtraction constants «o, /?o are fixed by the requirement 



M Q (s A ) = M^ oop 



-^M (s A ) = ^-Ml- Xoo *{s A ). (6.12) 

rl.ft- n.fi- 



ds a ds a 

The 1-loop amplitude enters here only to fix the subtraction constants at a well 
suited kinematical point, the Adler point. As mentioned, the corrections at this 
point are very small so that this procedure is reasonable. There is however a further 
implicit assumption in the method, namely the choice of the 7 = 1,2 amplitudes. 
Although this choice affects the total amplitude, say at the center of the Dalitz 
plot, rather substantially, the resulting uncertainty on c is not important. The 
contribution to c can then be written as a integral over the discontinuity in the s a 
channel: 

gC*0 = 1 f da'D(B')f§(s')[a + f3(s' - s A ) + jM ( a ')] 



AuD(s A )J (s'-s A ) 



Solving eq. ( 6.11| ) by iteration we obtain & Mq > = (7.9 + i0.5)GeV 4 , which has to 
be compared to the contribution of Mq~ 1o ° p , c ( M o,i-ioop) _ 5 g5G e v~ 4 . Thus 

5c = 2.05GeV~ 4 (6.14) 

which lies perfectly well within the range established above. However, the sign 
of the correction (positive) is now well understood. The reason is that current 
algebra underestimates the 1 = s-wave rescattering, hence yielding too small a 
discontinuity. 

We conclude that the actual value of 5c is higher than the one-loop value; and we 
will adopt c( M °) = (7.9 + i0.5)GeV~ 4 as the most likely number. It shifts the ampli- 
tude at the center to 1.71 -MO. 45 (of which 1.47 is in the subtraction polynomial) 
or by about 9 % []. 

3. Asymptotic behaviour of the 7T7r-phase shift 

In sect. 4 we mentioned an ambiguity to the solution of Khuri-Treiman equations 
due to the so far unspecified asymptotic behaviour of the 7T7r-phase shifts. Here we 
give a quantitative estimate of the effect on our results. We neglect the 1=1,2 



'We have taken the subtraction point SP = 0. 
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phase shifts and consider the approximate solution F = P + $10, where $10 is 
given explicitly in eq. ( |5.28 ). This approximation is actually very close to our final 



numerical solution, c.f. Table 3. Now we study two cases with distinct asymptotic 
behaviour of both, phase shift and amplitude: i) 5q(s) — ► and Q — 0, i.e. F(s) ~ s 2 
and ii) 5q(s) — > n and Q = const. 7^ 0, where we fine tune Q such that the amplitude 
has the improved asymptotic behaviour F ~ s. In the first case we use the phase 
shift as employed before, given in eq. (|B.2| ) and denoted here by 5q' Schcnk (s). In 
case ii) we need a model which guides the phase to its asymptotic value tt. We take 

frO, Schenk/ \ ^ 
% {s), s < S 0Q 

^ SchCn \s)(£ h , s 00 <s<Al 



71. 



s>A 2 h , 



with soo = (865MeV) 2 and A\ is the scale where the phase reaches its asymptotic 
value. We consider two choices h — 1/4 and 1/8 with corresponding scales A\u = 
(2.16GeV) 2 and Aj /8 = (4.87GeV) 2 respectively. 

We then calculate the difference 5F = — at the center of the physical 
region for various subtraction points SP. The maximal difference we obtain is 
5F = 0.015—2*0.020, or 1 % and 5 % of the full amplitude for real and imaginary part 
respectively. This in turn implies an error on the rate T of less than 3 % or w 4eV. 
The exercise shows that rather different assumptions on the asymptotic behaviour 
imply only modest changes of the amplitude in the low energy region. Here, taking 
three subtractions pays visibly. The same estimate for the twice subtracted case 
yields an uncertainty in the real part of the amplitude of already ~ 5% or m 17eV 
on the rate. 



4. Technical uncertainties 

We continue with some remarks on the structure and convergence of our iteration. 
As described in section 5, we distinguish two iterations, namely the one in powers 
of the weaker I = 2,1 rescatterings (^-iteration) and the iteration when solving 
the dispersion relation (m- iteration). In fact, the latter is an iteration in the 1 = 
rescatterings. We will discuss some exemplary numerical results and comment on the 
general trends. We will always use the subtraction at zero with the corresponding 
values of the polynomial. 

To illustrate the m-iteration, we take the succesive values of the contribution to 
pure 1 = problem as considered by Neveu and Scherk. Writing the amplitude U 
in the form (see eq. ( |5.17| )) 

00 

U k {s)=J2$mk(s) (6.15) 
1 

we find in the center of the Dalitz plot the following numbers: 

Re$i = 0.1828 Re <Z> 20 = -0.0631 Re $30 = 0.0051 Re $40 = 0.0001 (6.16) 
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This suggests that the iteration in the 1 = scattering phase converges nearly like 
a geometrical series with a coefficent of 0.3. We have found this behaviour over the 
whole Dalitz plot and for all values of k. The iteration is therefore terminated at 
m = 5. 

Next we consider the ^-iteration. As example, we take again the amplitudes U, 
V, W at the center. We obtain for the corrections to the polynomial: 

U = 0.1248 + i0.4786 no I = 1, 2 corrections 
U x = -0.0180 + i0.0638 one I = 1, 2 iteration 
U = 0.1145 + zO.5431 total result 
Vi = 0.0407 - i0.0617 one 1=1,2 corrections 
V = 0.0343 - i0.0649 total 

W 1 = -0.0038 - i0.0003 one 1=1,2 corrections 
W = -0.0039 + zO.0002 total 1=1,2 iteration 

Again we see that the convergence of the results is rather good; consecutive terms 
decrease by about 10 to 15 percent. 

To appreciate the degree to which the final result fulfills unitarity, we consider 
again the amplitude at the center of the Dalitz plot for various stages of the iteration. 
Including only the first iteration, # 10 , the value of the amplitude is A app = 1.57 + 
iO.5017. Unitarity is checked by inserting this (and the corresponding values for all 
s) value into the dispersion relation for the total amplitude which yields a new value, 
say A' app on the left hand side. Then the relative difference d = (A app — A' app )/\A' app \ 
is calculated. At the center of the Dalitz plot, we obtain for this deviation 0.0441 — 
z0.0301 which indicates that unitarity is fulfilled to about six percent. Taking now 
the full amplitude, we get at the center again 1.57 + z0.4128. The unitarity check 
now yields d = 0.0007 + i0.0002 for the deviation, a rather satisfactory result. We 
see that the iteration has improved the unitarity drastically, although the central 
value of the amplitude is not much changed. 

Another result of interest is the quality of the Neveu-Scherk $io iteration and 
the relative sizes of the various isospin amplitudes U, V and W. As seen from 
the values given before, P + $io corresponds rather accurately to the final value. 
This is accidental. If we just solve the complete Neveu-Scherk problem, e.g. only 
include 1 = rescatterings, the result would be U = 0.1248 + i0.4786 and $i 
would not be a good approximation. On the other hand, the complete calculation 
yields U = 0.1145 + i0.5431, V = 0.0343 - z0.0065, W = -0.0039 + z0.0002. This 
shows that the contribution from V is quite substantial, and that U alone does not 
saturate the amplitude at all. Nevertheless we note that \V/U\ is roughly about a 
third, which justifies our perturbative treatment of V (and W). 

Our calculation further depends on input parameters which are not well known. 
One is the cutoff of the dispersion integral. Changing it from 82m^ to 164m^., the 
rates vary at most by 1 eV and thus the variation is negligible. 

Also, the low energy constant L 3 |2j is not very precisely known. Varying L 3 
within the error bars we find the following ranges 
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r^+K-^o : 182.1 - 181. leV 
r^^^o^o : 250.9 — 253. 5eV 

a : -1.21 1.12 

b: 0.25-0.23 
c: 0.10-0.09 

for the one-loop value of the constant c. 

Apart from the changes in a, these variations are small. A more precise measure- 
ment of a could therefore restrict L 3 . 



6.2 Rates and Dalitz Plot Distribution 

Let us turn now to the values of the physical quantities. If we take the one-loop 
value for the constant c, the rate for 7] — > 7r + 7r~7r°, normalized by ml, is 

r = (180 ±40) eV. (6.17) 

The error given reflects our estimate of the uncertainty of 25 % of c. On the other 
hand, if the improved value for c is used, the rate becomes 

T = 209 eV. (6.18) 

We estimate the remaining error on this result to be about 20 eV. Since the one- 
loop value of c is most certainly too small, the result in eq. ( |6.18| ) should be more 
reliable than the value in eq. (|6.17|) The corresponding values for the decay into 



neutral pions are (252 ± 56)eV and 295eV. Thus, the rates are not substantially 
larger than the previous results, if the the one-loop value of c is taken. However, 
if its improved value is employed by taking into account the rescatterings through 
the Khuri-Treiman equations, the corrections are substantial, although the result 
is still considerably below the experimental one. The size of the corrections are 
qualitatively reasonable: since the one-loop corrections amount to about 50 % in 
amplitude, next order corrections might be expected at about 40 eV, roughly the 
value we found. 
The ratio 



Pjj— »7r 7r c) 7r () 

Ffj — >7r+7T~ 7T" (J 



(6.19) 



between the rates of the neutral and charged pion channel takes the values 

r = 1.40 ±0.03 (6.20) 

and 

r = 1.41 ±0.03 (6.21) 

for the two values of c, respectively. The values correspond to the the small subtrac- 
tion points SP considered; the error is obtained by varying the SP between —10 
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and 2 (see Table 4). In comparison, the value r = 1.43 was obtained in Ref. |J. 
The particle data group quotes 



r = 1.35 ±0.05 "our fit" 

r = 1.27 ±0.14 "our average", (6.22) 

favoring smaller values for r, in particular almost excluding the current algebra 
result. Recently, r has been remeasured in a direct measurement P7| 



r = 1.44 ±0.09 ±0.01 (6.23) 

with smaller errors than previously and which is in very good agreement with our 
result and the one-loop value but somewhat below the current algebra prediction of 
r = 1.51. 

The distribution over the Dalitz plot is conventionally described by the two vari- 
ables x and y 

3 

V = 2m Q ^ mri ~ ~ Sa) ~ 1 < ^ 6 ' 25 ' ) 
Q v = m r] — 1m 1 s+—m^. (6.26) 

Since Q v is small, it is important to keep the masses of the charged and neutral 
pions different. The mass difference is generated by electromagnetic corrections. 
However, as discussed in the first section, the other electromagnetic contributions 
to the decay amplitude are small and can be neglected. 

The Dalitz plot distribution \A(x, y)\ 2 can be parametrized by a second order 
polynomial. For the decay into charged and neutral pions we may write, respectively, 



\A + -°(x,y)\ 2 = N + -° (l + ay + by 2 + cx 2 
\A 000 (x,y)\ 2 = N 000 (l + g{x 2 + y 2 )). (6.27) 

Of course, the neutral pion mass must be used for x and y in the second expression. 
The various determinations of the slopes a, b, c and g are as follows: 
Gasser and Leutwyler §: a = -1.33, b = 0.42, c = 0.08 

Layteretal. (exp)[p|: a = -1.08 ± 0.014, b = 0.034 ± 0.027, c = 0.046 ± 0.031 
Gormleyet al. (exp)|||: a = -1.17 ± 0.02, b = 0.21 ± 0.03, c = 0.06 ± 0.04 
Amsleret al. (exp)[0: a = -0.94 ± 0.15, b = 0.11 ± 0.27 
Aide et al. (exp)gl|: g = -0.044 ± 0.046 

This work: a = —1.16, b = 0.24, c = 0.09, g = —0.028 (one-loop value for c) 
a = —1.16, b = 0.26, c = 0.10, g = —0.014 (improved value for c). 
We see that our numbers are very near to the experimental values, in fact closer 
than the one-loop results. The values of a for negative SP are somewhat too large, 
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while the result for SP = correspond nicely to the result of Gormley et.al. Unfor- 
tunately, the two most accurate experiments are mutually inconsistent and thus pre- 
clude a definite statement; measuring the slopes with larger accuracy would remove 
the uncertainties inherent in the choice of the subtraction points. The quadratic 
slope parameter of the neutral decay mode, g, is experimentally compatible with 
zero. We observe that the slopes (with the exception of g) do not depend strongly 
on the value of c. 

The shape of the amplitude, however, changes compared to the one-loop ampli- 
tude. In Fig. 3 we plot the real part of our numerical solution (with the onel-loop 
value of c) along the line s a = s c , together with the current algebra prediction as 
well as the chiral perturbation theory one-loop amplitude. At small s, all amplitudes 
are close together which reflects the fact that we have subtracted at s a = 0. The 
cusp generated by the two-pion threshold is more pronounced in the solution taking 
into account the corresponding final state interactions to all orders. However, the 
requirement of unitarity bends the amplitude down more strongly in the physical 
region, leading to a value of the amplitude at the center of the Dalitz plot close to 
the value obtained to one-loop ChPT. At the same time, the slope of the amplitude 
is reduced substantially which yields a smaller values of the linear slope parameter 
of the Dalitz plot distribution. The imaginary part is seen to be enhanced over the 
whole physical region, c.f. Fig. 4. This was expected, as the one-loop amplitude 
gives only the leading term to the imaginary part. The figures for the improved c 
have a similar shape, but with different normalization corresponding to the larger 
rate discussed above. 



6.3 Two Subtractions 

We turn briefly to the twice subtracted dispersion relation. Starting from the current 
algebra polynomial, we obtain the results in Table 6, taking the subtraction points 
in Table 5. We have used the same scattering phase as in the three times subtracted 
case, and set the polynomial Q equal to zero. This time, the rates grow fast if 
we change the subtraction points to large negative values. On the other hand, 
this increase is accompagnied by an unreasonable change of the slopes a and b 
as well as unacceptable values of r. We expect similar changes, if we vary the 
asymptotic behaviour of the amplitude and the phase. We conclude that artificially 
enhancing the rates with large negative subtraction points does not yield correct 
results and that the higher order terms in the subtraction polynomial are essential. 
An alternative approach would be to use the one-loop subtraction data, but maintain 
a twice subtracted relation. We expect in this case larger uncertainties from the 
assumption of elastic unitarity as well as from the asymptotic ambiguities. 



Roiesnel and Truong |32] were the first to invoke unitary corrections to enhance 
the decay rate. These authors obtained satisfactory values for the rates, using a once 
subtracted dispersion relation; furthermore they only considered 1 = rescatterings, 
but only for one part of the amplitude. As a result, their discontinuity is not given by 
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the usual form, but larger by the factor 9/5. Since the sign of the 1 = 2 contribution 
is negative, this factor in fact gets enhanced because the I = 0, 2 pieces suffer an 
artificial cancellation. If this factor is omitted, the result is not as dramatic and close 
to our results for SP = 0. As noted, the dispersion relation with one subtraction is 
subject to large ambiguities and thus the result will be beset with high uncertainties. 



7 Conclusions 

The decay amplitude of rj — > 3tc is proportional to the quark mass difference (rrid — 
m u ) and thus provides one of the ingredients to determine the important ratio ^ 2t . 
Including also higher orders in chiral perturbation theory, the decay rate can indeed 
be written as 

r = {Qdt/QYT (7.1) 

where 

_ m d -m u m d + m u 
m s — rh m s + rh 

with rh = + m u ). Qdt is the value of Q if Dashen's theorem for the electro- 
magnetic kaon mass difference is used (see eq. ( p.. 9ft ) and f the corresponding width 
as given in the previous chapter. 

Using the value Qdt (or the corresponding and long established values for the 
quark mass ratios), the one- loop ChPT prediction for the rate given by Gasser and 
Leutwyler || has been considerably below the experimental number. Theoretically, 
the rate can be increased by lowering the ratio m u /md or by including fur- 
ther corrections in f. In particular, it has long been suspected that the unitary 
corrections may be sufficient to enhance the rate sufficiently 



The existing framework embodied in the so called Khuri-Treiman equations [24 



allows to take into account rescatterings of two pions which are thought to dominate. 



Some time ago, Roisnel and Truong JJIJ claimed that in fact these corrections solve 
the problem, but as explained in section 6 (see also 0), we believe that their results 
are an overestimate. 

In this paper we have calculated the unitary corrections using the Khuri-Treiman 
equations (which include the I = 0, 2 pion pion rescatterings), complemented with 
the 1 = 1 interactions. To specify these dispersion relations completely, that is 
to define the subtraction method, we have used an iterative procedure for solv- 
ing it which accounts in lowest order for the one-loop ChPT results of Gasser and 
Leutwyler. This forced us to a three times subtracted dispersion relation with rather 
fast convergence and little dependence on the unknown high energy scattering phase. 
In this approach, the subtraction polynomial is quadratic in the invariant momenta, 
rather than the linear form of current algebra. The major source of errors lies in the 
choice of the subtraction points and, particularly, in the uncertainties in the subtrac- 
tion constants. The constants connected to the quadratic terms vanish at tree level 
and only start at the one-loop level; thus their one-loop value is rather uncertain. 
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In section 6 we have described how this problem can be overcome (following recent 
work of Anisovich and Leutwyler f22|), and have given a reasonable estimate of the 
constant c. We have chosen to give the results corresponding to both their one-loop 
value and to the improved one. Using everywhere the new value of f„ |J, the rate 
corresponding to the the one-loop value is 

f = (180±40)eV, (7.3) 

while the improved number of c yields 

f = 209eV (7.4) 

and where we estimate the remaining errors to be about 20eV. We view the second 
result as the more reliable one. Although the corrections are quite large, they do 
not suffice to explain the experimental value of (281 ± 28)eV. 
From eq. ( |7.4| ) we obtain 

Q = 22.4 ±0.9. (7.5) 

while the rate of 180 eV in ( |7.3| ) would imply Q = 21.6±1.3. In contrast, Qdt = 24.1. 

These lower values of Q can obviously be accounted for if the ratio ^ is lowered, 
and another mass ratio, say md+mu is kept fixed. In this case, — is reduced from 0.57 



to 0.49 ± 0.04 or to 0.52 ± 0.03 depending on whether eq. (|7.3|) or eq. ( |7.4|) is used. 



Thus, while it may be somewhat smaller than previously thought, the up quark 
mass does not vanish FI. We note that the somewhat lower value of — in l20| 

I — I rriri II U 



was 



obtained with a changed value of md " . The constant Q can also be determined 
from electromagnetic corrections to the meson masses. As already noted in |20|| , 



these lower values of Q correspond roughly to the one obtained from the kaon mass 
difference, if electromagnetic corrections are positive and large as found in (TH . 



We have also determined the slopes in the Dalitz plot distributions; the results are 
given in the last section. We find that the ratio r between the rates of the neutral 
and charged pion mode remains roughly at 1.4 which is also favored experimentally 



37] . On the other hand, the slopes change from their one-loop ChPT values; our 
number for a is higher than the previous result and nearer to the experimental value. 

The experimental situation is rather unsatisfactory. The decay rates are normal- 
ized with respect to the decay with T(rj — > 77), where the results from two-photon 
production disagree with those from Primakoff production. The two-photon mea- 
surements seem more reliable; however, in order to resolve the issue completely a 
reanalysis of the Primakoff data would be necessary. The Dalitz plot distribution of 
i] — > 7T + 7r~7r° has been measured with rather high accuracy |59|, [SB]]. However, the 
assumptions made by these authors are not compatible and do not allow comparison 
ot the numerical values. A recent experiment [B7|] has still too large error bars in 

10 our analysis is based on the assumption of the validity of standard chiral perturbation theory. 
In contrast, a treatment along the lines of Ref. jl4j might lead to different results; however the 
larger quark masses in that scheme exclude a zero mass automatically 
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order to be conclusive. As our results show, there is a rather strong correlation 
between rates and slope parameters and a more acurate measurement of the latter 
would fix the rates, and thereby the up quark mass better. 
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Appendix A 

In this Appendix we give the three projections of the subtraction polynomial P. 
Since P has no cuts, the integrals as defined in section 3 are trivial and may be 
performed explicitly. 

The second order subtraction polynomial is of the form 

P(s a , s b , s c ) =a + (3s a + js 2 a + 5(s b - s c f (A.l) 

and we need the projections for n — 0, .., 3. The corresponding results are 



1 



1 

3s - s a 



3s -s a \ 2 K 2 (s a ) 
2 J 3 
3s - s a 



3s - s a 



+ K 2 (s a ) 



We now insert P into eq. ( |3.18|) , reexpress s c = 3sq — s a — s& and obtain 



P(s a ) = a + (3s a + js 2 a 

+ 5{4~4 - A- b {3s - s a ) + (3s - s a ) 2 ) 

= a + f3s a + 1 s 2 a + 6^P±. 



(A.2) 



(A.3) 



For P(s a ) we get from eq. (|3.19|) 



P(s a ) = a + Psb + 'ysl 



5(s 2 - 2- b {3s - 2s a ) + (3s - 2s a ) 2 ) 
a + /3 3g ° - 5s a (3s - 2s a ) 



+ (7 + 5) 



/ 3gQ - Sg 

V 2 



(A.4) 
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p ( s -) = - 7^7: A 2{as- b + Psl + 1S 3 b 



Turning finally to P its insertion into eq. ( |3.30| ) yields 

3 

+ 5(4-24{3s -2s a ) + Sb{3so-2s a ) 2 } 

-(3s - s a ){a + f^ + 'jsl (A.5) 
+ 6(4 - 2s- b (3s - 2s a ) + (3s - 2s a ) 2 )}] 



~\J3 + { 7 -S)(38 -s a )]. 



Appendix B 



In this Appendix we discuss the parametrization of the different pion-pion phase 
shifts used in our numerical computations and disply the resulting Omnes function. 

A careful analysis of the different restrictions on the phase shifts is carried out 
in Ref. There are two main points to be respected, namely first the threshold 
behaviour of Re// 

Re/fa) = ( S -^Y (a\ + 6f(^=i) + ...) (B.l) 

where we used the notations of section 3 and where af denote the different scattering 
lengths, b\ the corresponding slope parameters. Second, one has to implement that 
the phase shifts pass through | at some experimentally known values s = sf of the 
energy. A simple parametrization respecting those conditions is given by |^1 



tan ^ (s) = ^i_i.(l_i ) ..|_i (B.2) 

where the threshold expansion is reproduced with 

4 

4 

As the experimental data on pion-pion scattering near threshold are rather poor one 



b'^bj-aj-^ + ia^o. (B.3) 



uses the results of ChPT for the scattering length [13, 42 1 



a° = 0.20, a 2 Q = -0.042, a\ = 0.037, (B.4) 



b° = 0.24, b 2 Q = -0.075. 



The remaining data are extracted from experiment 

b\ = 0.005, 

c o = , c 2 = 0, c} = 0, (B.5) 
s° = 38.45, s 2 = -24.11, s{ = 30.39 
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where sf is given in units of m 2 as usual. These values correpond to the 'best fits 
in 



41fl . We remark that the constraints of the Roy equation around s = 4 are taken 
into account in this parametrization. 

We already noted that all the numerical integrals will be cut off. In order to avoid 
discontinuities in the numerical integrations, caused by the step-function induced 
from the cut-off, we bring the phases smoothly down to zero multiplying them with 
the exponential suppression factor 

e~ k < s - s ^ 4 (B.6) 

ensuring differentiability at Sq. For k we chose the value k = 10~ 5 . 

We do not give explicitely the numerical result for the Omnes function. We just 



note that instead of working directly with the definition in eq. ( |5.10| ) we prefer to 
go over to the once subtracted form 

. r°° *-> s ° (3,) 

D(z) = e A^). (B.7) 

This is possible as we are dealing only with fractions of such functions in which the 
extra contribution from the subtraction point at z = cancels. The advantage of 
subtracting here is the better numerical convergence of the integrals. 



Appendix C 

In this Appendix we analyse possible consequences of the singularities occuring 
in the projected amplitudes as displayed in eqs. ( |3.38| ) and ( |3.39| ). We show in 



particular that for physical values of the energy, s + ie, the amplitude resulting from 
the dispersion integrals is finite. 

Within the iterative procedure we have to perform dispersion integrals over pro- 
jected functions having poles at s = (m — l) 2 to obtain certain contributions to the 
amplitudes in question. As long as the value of the external variable s is not around 
s = (m — l) 2 there are no problems with the integrations. We have thus to analyse 
only the case s ~ (m — l) 2 where the two types of integrals are of the form 

,(m-l) 2 -S+g I I 

const. ■ / ds 



(m _i)2_ 5 _ 3 s'- (m- I) 2 + 5 -is (( m _ i)2_ s /) 5 + n 

+finite. (C.l) 

Above we introduced the two positive constants g > 5 > dealing thus with the case 
s = (m — l) 2 — 5 — > (m — l) 2 from below. The other case is then treated similarly. 
g should be small in the sense that we may replace the full integrand by its above 
reduction to the two rational functions, 'finite' denotes the obviously finite rest and 
n the two cases coming from eq. ( |3.38|) with n = and from eq. ( |3.39| with n = 1. 
After a change of variable s' — > x = s' — (m — l) 2 + 5 we obtain 

r+9 1 1 

I n = / dx — (C.2) 

J -9 x-ie (- x + S)h+ n 
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discarding all unimportant contributions to eq. ( p.l| ). Using the Sokhotsky-Plemelj 
formula we obtain the principal value integral 

l r+9 1 1 

I n = m^— + P dx j—. C.3 

As the square-root has a cut we have to distinguish the cases — x + 5 > and 
-x + 5 < 0. 

For —x + 5 > we obtain 

dx- — i — = [ dz^- \- (C.4) 

x (- x + S)h+ n 52 +n J z 2 -l z 2n v ; 



making the substitution z = J x ^ 5 . For n = this leads to the contribution 

-log ^4 , (C.5) 



52 9+1 + l 



for n = 1 to 



+ log ^J- 1 - 2 • lim J \ ( ( '.(i ) 



5§ 



where we properly distinguished the two cases \z\ < 1 and \z\ > 1 in the range of 
integration and where the principal value prescription leads to the cancelation of 
the infinities at x = 0. 

For —x + 5 < we have to continue analytically. With \/—x + S = iy/ x — 8 we 
find 

x— <5 



making the substitution z = \J . For n = this leads to the contribution 
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for n = 1 to 



Cretan, fc^, (C.f 



2i 5 g-d , 6 

+ arctan \ — hm W > . (^-9) 



5^ I V 9 ~ ^ V 8 x^l \ x — 5 



Here no further distinctions have to be made. 
We collect all the contributions for n = 



1 f /— — 1 / ~5 

I = — I in — log 2i arctan \ — - — \ (CIO) 
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and for n = 1 



h = (Cll) 




— 2i arctan y - - - + 2 • lim 

In the limit 5 — > a Taylor expansion shows that the log-terms vanish whereas the 
arctan-terms cancel the ^-contributions such that I remains finite. 

In the case n — 1 there remains the divergent contribution ~ Resubstitu- 
tion of x = s — (m — l) 2 + 5 yields a square-root singularity at s = (m — l) 2 . As 
this type of behaviour occurs in the function W(s a ) only which is always multiplied 
with Sb — s c = 2cos6 ha ■ K(s a ) we see that the square-root singularity is lifted by 
the Kacser function K. The case s = (m — l) 2 + 5 leads to the same results in the 
limit 8 — > 0. The resulting amplitudes are thus finite at s = (m + l) 2 + is although 
discontinuities may occur at this point. 

At the unphysical boundary s = (m + l) 2 — is the amplitudes badly diverge as 
the arctan contribution does not cancel anymore the 27r-terms as above. 



Appendix D 

Here we give the explicit form of the functions RS 1 ^ occuring in the determination 
of the subtraction polynomial for finite subtraction points, c.f. section 4, eq. (|4.31|) . 
The derivation employs the identity 

ReAj(s) = (s - Si)(s - Sj)y(s,Si,Sj) 

+ S - [ Si (y(0, Si , Sj ) + y(0, 0, s< )) + (i <-> j)} 
+s i s j y(0,s i ,s j ), (D.l) 

A 7 (s) = - / ds'-Z-, (D.2) 

V ' 71 Jiml S' 2 { S >-S-ie) V ' 

1 ,oc X/,ChPT/ A 

y I (s,s i ,s j ) = -P ds'- y 4r 7- ( D -3) 

y V ' ' 3) 7T hml {S' - S){S' - Si){s' - Sj) V ' 

P denotes the principal value and J / > ChPT is the 7T7r-phase with isospin / to leading 
order in ChPT. For / = 0, 2 we thus find 

RV\s; Sl ,s 2 ,s 3 ;s A ) = ^—11 £ || [* (V (0, Si , 8j ) + y J (0, 0, *)) + (z <- j 



where 



and 



/(OjSijSi)} 
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1,2,3 



+y I (s 1 ,s 2 ,s 3 )- ( s - s i)( s ~ s j)( s A~ S k ). (D.4) 



i<k,k^i,j 



The corresponding expression for / = 1 is 
R {1 \s a ,s b ,s c ;w 1 ,w 2 ) = 



x 



^ 2 ^ ^2 / \ 2 

vn (V (0, W!, w 2 ) + y(0, 0, Wl j) + (1 ~ 2) 



+ 2~( Sa ~~ s )w 1 w 2 y 1 {0,w 1 ,w 2 ). 



(D.5) 
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Tables 



SP 




<s 2 


«3 


-10 


-9.6 


-9.8 


-10.0 


-6 


-5.6 


-5.8 


-6.0 


-4 


-3.5 


-3.7 


-3.9 


-2 


-1.7 


-1.9 


-2.1 





5* 1(T 5 


1.2* 10 -3 


5 * 10~ 3 


1 


1.2 


1.33 


1.4 


2 


1.7 


1.9 


2.1 


3 


2.6 


2.7 


2.8 



Table 1: Values of the subtraction points SP used in the three times 
subtracted dispersion relation in units of m^. The subtraction points 
are the same for U and for V, for W we take the first two values on each 
row (see also text). 



SP 


a 


P 


7 


5 


-10 


-1.38 


19.73 


-7.87 


3.74 


-6 


-1.28 


20.52 


-5.45 


4.00 


-4 


-1.25 


20.97 


-3.45 


4.13 


-2 


-1.27 


21.39 


-1.90 


4.22 





-1.28 


21.81 


4.09 


4.19 


1 


-1.34 


21.99 


11.24 


3.88 


2 


-1.38 


21.90 


17.46 


3.44 


3 


-1.47 


21.08 


36.67 


1.81 


error bar 


±0.14 


Tl.52 


±3.18 


=Fl.08 



Table 2: Coefficients of the subtraction polynomial for the sets of sub- 
traction points specified in Table 1. The last line gives the error bar due 
to the uncertainty in the determination of L 3 . All values given here and 
in the subsequent tables refer to the one-loop value of the constant c 
discussed in section 6. 
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SP 


P 


ReGL 


A 


Act 


-10 


0.86 


0.64 


1.62 + i0.44 


1.63 + «0.41 


-6 


1.10 


0.40 


1.70 + i0.49 


1.69 + ^0.44 


-4 


1.21 


0.29 


1.69 + i0.50 


1.67 + i0.44 


-2 


1.26 


0.24 


1.60 + ^0.49 


1.59 + ^0.42 





1.40 


0.10 


1.57 + i0.50 


1.57 + «0.41 


1 


1.46 


0.11 


1.51+^0.40 


1.50 + ^0.39 


2 


1.49 


0.01 


1.42 + «0.49 


1.46 + ^0.38 


3 


1.58 


-0.08 


1.32 + ^0.48 


1.41+^0.36 



Table 3: The polynomial approximation P and the real part of the dis- 
persive part of the one- loop result, ReGL, in the center of the Dalitz 
plot. The sum is fixed to give 1.5, the value of the one-loop amplitude. 
Also the result of the first approximation (including only $ 10 ) and the 
total amplitude are given 



SP 


a 


b 


c 


9 






r 


-10 


-1.27 


0.32 


0.07 


-0.023 


201eV 


278eV 


1.38 


-6 


-1.22 


0.29 


0.08 


-0.020 


213eV 


297eV 


1.40 


-4 


-1.21 


0.30 


0.08 


-0.017 


211eV 


295eV 


1.40 


_2 


-1.21 


0.30 


0.08 


-0.020 


189eV 


2Q3eV 


1.40 





-1.16 


0.24 


0.09 


-0.028 


180eV 


253eV 


1.40 


1 


-1.04 


0.13 


0.10 


-0.046 


162eV 


2UeV 


1.44 


2 


-0.97 


0.05 


0.09 


-0.082 


U8eV 


2UeV 


1.45 


3 


-0.78 


-0.06 


0.08 


-0.107 


128eV 


189eV 


1.48 



Table 4: Values of the physical quantities for different subtraction points 
(three times subtracted dispersion relations). 



SP 


Sl 


S2 


-10 


-9.8 


-10.0 


-6 


-5.6 


-6.0 


-4 


-3.7 


-3.9 


-2 


-1.9 


-2.1 





1.3* 10" 3 


5* 10~ 3 


1 


1.2 


1.33 


2 


1.9 


2.1 


3 


2.7 


2.8 



Table 5: Values of the subtraction points SP in units of ml for the twice 
subtracted dispersion relation 
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or 


a 


u 




c 




1 000 


r 


-10 


-0.81 


0.12 


0.04 


410eV 


Q2AeV 


1.52 


-6 


-0.72 


0.10 


0.06 


277e\/ 


423eV 


1.53 


-4 


-0.77 


0.14 


0.05 


222eV 


335eV 


1.51 


-2 


-0.90 


0.20 


0.06 


182eV 


2Q8eV 


1.47 





-1.14 


0.37 


0.03 


U5eV 


203eV 


1.40 


1 


-0.96 


0.28 


0.04 


120eV 


173eV 


1.44 


2 


-0.78 


0.20 


0.05 


95eV 


UOeV 


1.47 


3 


-0.71 


-0.13 


0.02 


78eV 


117 'eV 


1.50 



Table 6: Values of various physical observables for different subtraction 
points and two subtractions. 





5c/ 'c = +25% 


Sc/c = -25% 


5d/d = +25% 


0.15 + i0.04 


-0.13-i0.03 


5d/d = -25% 


0.18 + i0.05 


-0.09-i0.02 



Table 7: Change of the reduced decay amplitude A at the center of the 
physical decay region for a relative error of 25 % on constants c, d. 



Figure captions 

Fig. la Integration contour for case i) 

Fig. lb Integration contour for case ii) 

Fig. lc Integration contour for case iii) 

Fig. Id Integration contour for case iv) 

Fig. 2 Lattice of points for the numerical determination of the different functions 

Fig. 3 Real part of the decay amplitude for s = u (in units of m^) for SP = 0. 
The full line is the numerical solution of the coupled integral equations as 
described in the text. Also shown is the current algebra prediction (dashed) 
and the chiral perturbation theory one- loop result (dotted). The physical 
region lies between the two vertical lines. 

Fig. 4 Imaginary part of the decay amplitude for s = u (in units of m^) for 
SP = 0. The full line is the result of the numerical iteration and the dotted 
one the chiral perturbation theory one-loop approximation. 
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Fig. lc 



Fig. Id 



/ 


1 








4 m+1 



Fig. 2 




Figure 3 



physical region 



I I I I I I I I I I U I I I U I I L 

2 4 6 8 

s 



Figure 4 



